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ABSTRACT 

Aims. Our aim is to characterize the structure of protostellar envelopes on an individual basis and to correctly identify the embedded 
YSO population of L 1688. 

Methods. Spectral maps of the HCO+ / =4-3 and C'*0 7 =3-2 lines using the HARP-B array on the James Clerk Maxwell Telescope 
and SCUBA 850 fim dust maps are obtained of all sources in the L 1688 region with infrared spectral slopes consistent with, or close 
to, that of embedded YSOs. Selected 350 yum maps obtained with the Caltech Submillimeter Observatory are presented as well. The 
properties, extent and variation of dense gas, column density and dust up to 1' (~ 7,500 AU) are probed at 15" resolution. Using the 
spatial variation of the gas and dust, together with the intensity of the HCO^ J =4-3 line, we are able to accurately identify the truly 
embedded YSOs and determine their properties. 

Results. The protostellar envelopes range from 0.05 to 0.5 M© in mass. The concentration of HCO^ emission (~ 0.5 to 0.9) is gener- 
ally higher than that of the dust concentration. Combined with absolute intensities, HCO^ proves to be a better tracer of protostellar 
envelopes than dust, which can contain disk and cloud contributions. Our total sample of 45 sources, including all previously classified 
Class I sources, several flat-spectrum sources and some known disks, was re-classified using the molecular emission. Of these, only 
17 sources are definitely embedded YSOs. Four of these embedded YSOs have little (0.1-0.2 Mq) envelope material remaining and 
are likely at the interesting transitional stage from embedded YSO to T Tauri star About half of the flat-spectrum sources are found 
to be embedded YSOs and about half are disks. 

Conclusions. The presented classification method is successful in separating embedded YSOs from edge-on disks and confused 
sources. The total embedded population of the Ophiuchus L 1688 cloud is found almost exclusively in Oph-A, Oph-B2 and the 
Ophiuchus ridge with only three embedded YSOs not related to these regions. The detailed characterization presented will be neces- 
sary to interpret deep interferometric ALMA and future Herschel observations. 



• ^ ' 1. Introduction Tauri stars with gas-rich circumstellar disks and Class III the 

, „ ,. X , • , pre-main sequence stars surrounded by tenuous or debris disks. 

^ Low-mass young stellar objects (JSOs) have traditionally been photometry introduced the 'flat-spectrum' (FS) 

C3 ; classified using thei r observed in frared (IR) spectral slope ffjR, ^^^^^^^ ^ Greene et al.l fT994). with IR spectral slope close to 

from 2 to -20 ^m ( |Lada & Wilk ing 1984 Adams et al. 1987) q ^^-^^ represent a stage intermediate between Class I 

or their bolometric temperature, Tbo, ( Myers & Ladd| 09931) . „ ^^^^^^^^ classification and physical characterization 

Together wi t h the subseque nt discovery of^ Class stage ySOs is important for constraining the time scales of each of 

(| AndreetdJ U 993|) this led i GreeneetalJ to identify 5 ^^^^^ ^J^^^ determining the processes through which an 

object transitions from one phase to the next. In this study, we 

- Class , no az-u, high Lsub-mm/iboi f^'^'^' °" embedded YSO population and their transition to 

- Class I , ^2-14 > 0.3, Tboi < 650 K the T Tauri phase. 

- Flat Spectrum, -0.3 < 02-14 < 0.3, Tboi ~ 400-800 K 

- Class II -2 < a2-H < -0.3, 650< < 2800 K .^j^^ ^ Ophiuchus molecular clouds, part of the Gould 

- Class III , a2-i4 < -2, Tboi > 2800 K g^j^^ ^^^^ nearest star-forming regions and contain 
Each class is fliought to be represent a different category, and "^^"y Class I and II sources. Consisting of two main clouds, 
probably evolutionary stage, of YSOs. Class sources are ^ 1688 and L 1689, the star formation history and protostellar 
the earUest, deeply embedded YSO stage; Class I sources are population of fliese regions have been studied extensively, 
thought to be more evolved embedded YSOs, Class II the T Although the distance to Ophiuchus has long been debated 

(Knude & H0g 1998), recent work constrains it to 120+4 pc for 

Send offprint requests to: T. A. wan Kempen L 1688 ( Loinard et aLl l2008l) . The large-scale structure of the 

e-mail: tvankempen(acf a .harvard, edu Ophiuchus clouds at millimeter wavelengths was first mapped 
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by iLoreiil (119891) using the '^CO molecular line emission 
with 2.4' resolution. It was found that much of the cloud is 
in filamentary structures, but also that a diffuse foreground 
layer is present in Ophiuchus, resulting in a higher average 
extinction toward YSOs than in other clouds such as Taurus 
jPickman & Herbs^l 19901) . Subsequent continuum observations 
at millimeter (mm) wavelengths mapped most of the large-scale 
structure present in detail, distinguishing the Oph-A through 
Oph-F regions within L 1688 (Mezgeretal. 1992, Motte et al. 
19981 iJohnstone etsU l200a .Stanke et al., ,2006. .Young et al.i 



2006) 



The YSO populati on was fi rst identified b y lEliasI (Il978h. 
IWilking & Ladal (11983 ). Wilkin^etaD (Il989h . IComeron et all 
( 19931) aiid" Greene et al.i (il994i) using IR observations. With 
the arrival of (sub-)mm telescopes, VLA 1623 in L 1688 was 
identified as the first deeply e mbedded YSO dWoottenl [19891 
iLoren et"aL 1 119901 lAndreet al.l[r993i) . In more recent years, a 
large population of Class I and Class II sources has been found 
based on their IR spectral slopes, using space-based obser- 
vatories such as the Infrared Spac e Observatory (ISO) (e.g. 
iLiseau et al.l 11999 ', Bontemps et al.j 1200 lb and ground-based 
IR (e.g. iBarsony et al. 199 7; 2005h . Most embedded sources 
in L 1688 are clustered around the filaments of the Oph-A, 
Oph-B2, Oph-E and Oph-F regions, while the Oph-C region 
only shows a single em bedded source an d Oph-D does not have 
any embedded YSO ( Motte et al.lll998l) . In Oph-E and Oph-F 
most sources are lined up along a relatively small filament of 
material: for the purpose of this paper, we will adopt the name 
'Ophiuchus ridge' for this region. 



With the launch of the Spitzer Space Telescope, the 
Ophiuchus cloud was included in the guaranteed time (GTO) 
and t h e 'cores to d i sks' ( c2d) Legacy program (Eva ns et"al1 
l2003h . iPadgett et al] (l2008h report on the results at 24, 70 and 
160 fim using the MIPS instrument, revealing the emission 
from the large- s cale s tructure at mid and far-IR wavelengths. 
|j0rgensen et al.l (|2008") compared the results from the c2d pro- 
gram with the COM PLETE 850 yum SCUBA sub-millimeter dust 
mapping from IJohnstone et al.l ([2000) and Rid ge et al.l (i2006i) to 
determine the association of YSOs with dense cores. 

The stellar ages of the Class II and III sources in Ophiuchus 
were found to be 0.1-1 Myr based on stellar spectroscopy com- 
pared with evolution ary tracks, indicating a relatively young age 
for th e total cloud (iGreene & MeveJ I1995L iLuhman & Riekd 
Il999l) . The Star-Formation Efficiency (SEE) was recently 
calculated with Spitzer and SCUBA photometry to be of the 
order of 13% within the c ores and 4% in t he clo ud as found 
by .Evans et al.1 ( l2009l) a nd lJ0rgensen et all (l2()08l) . lower than 
previous determinations (I Wilkin g & Lada|[r983h . 



The relative timescales of the different phases are deter- 
mined by the nu mber of objects in each class of YSOs (e.g. 
lEvans et aLll2009l) . Recent high resolution ground-based (near)- 
IR imaging show that some of the Class I sources in Ophiuchus 
are physically different from an embedded YSO, confusing 
these timescales determinations. For example, the Class I source 
CRBR 2422.8-3423 w as found to be an extincted edge-on disk 
from near-IR imaging dBrandner et al.ll2000l iPontoppidan et al.l 
2005h. The source OphE MM3, classified as a starless core by 



Motte et al.1 (11998*). was also shown to be a edge-on disk in the 



same study (fBrandner et al. 2000). The source IRS 46 has no 
associated protostellar envelope and was re-classified b ased on 
Spitzer and sub-mm data as a disk (iLahuis et al.ll2006h . Much 



of the reddening seen in the IR originates from the nearby 
envelope associated with IRS 44. 

Foreground material can also heavily influence the iden- 
tification and sub sequent analysis of embedded sources (e.g., 
iLuhman & Riekei ri999). An excellent example is provided by 
the Class I source Elias 29 in L 1688, which has two fore- 
ground layers in a ddition to the ridge o f material in which the 
YSO is embedded dBoogert et al.ll200^ . Only a combination of 
molecular line emission at sub-millimeter and IR spectroscopy 
could constrain the protoste llar envelope as wel l as the imme- 
diate environment around it (Boogert et al. 2000; 2002). Indeed, 
IR spectroscopy can be used as a complementary diagnostic and 
spectra of many of the YSOs in Ophiuchus hav e been taken, us- 
ing ISO, Spitzer or groud-based telescopes (e.g. 'Ale xander et aTl 
2003, Pontoppidan et al. 2003, Boogert et al. 2008). Ice absorp- 
tion features such as the 3 jum H2O and 15.2 jjm CO2 bands 
are usually associated with embedded sources whereas silicate 
emission at 10 and 20 fj.m is characteristic of Class II sources, 
but foreground absorption and edge-on disks can confuse this 
classification (Boogert et al. 2002, Pontoppidan et al.ll2005h . 

In recent years, several detailed modelling efforts have been 
carried out to study the relations between the observed spec- 
tral energy distribu tion (SEP) and the physical structure of em- 
bedded YSOs (e.g. J0rgensen et al. 2002, Whitney et al. 20031 
Schoieret al. 200 4, Young et al. 2004, Robitaille et al. 200^ 
2007] iCrapsi et al. . 2008). Whit ney et al.1 (l2003bt) show that it 
is possible for embedded YSOs with a face-on projection to be 
classified as Class II. Due to their orientation, these sources are 
viewed straight down the outflow cone, directly onto the central 
star and disk system. Crapsi e t al. (2008) show that a significant 
fraction of the Class I sources may be edge-on flaring disks, that 
have already lost their protostellar envelope. The spectral slope 
is much steeper than expected due to their structure. 

A prime characteristic and component of embedded YSOs is 
the presence of dense centrally condensed envelopes. While dust 
maps at sub-millimeter wavelengths have be come very popular 
to trace the early stages of star formation (e.g. Motte et al. 199 
Shirley etal. 2000, Johnstone et al. 2000, Stanke et al. 200( 
the continuum emission at these wavelengths is dominated by 
the cold outer envelope and cloud material, with disks starting 
to contribute as the envelope disperses (e.g. iHogerheiide et al.l 
119971 iLooney etalj|2000l lYoung et al.|[2003h . Single-dish dust 
continuum data by themselves are not able to distinguish be- 
tween dense cores and envelope or foreground material, nor 
quantify any disk contributions. However, the dense gas (-lO*" 
cm"^) located in the inner regions of protostellar envelopes is 
uniquely probed by molecular lines with high critical densities at 
sub-millimeter wavelengths. Observations of deeply embedded 
Class YSOs have indee d revealed strong sub-millimeter Unes 
of various molecules (e.g.'Blake et al."l994';'l995L ISchoier et al.l 
2002, J0rgensen 2004, Maret et al. 2004; 2005), but only a few 
studies hav e been carried out on mo re evolved Class I embedded 
YSOs (e.g. IHogerheiide et al.|[T997h . 

A good high density tracer is the HCO^ molecule, for which 
the J =4-3 line both has a high critical density of >10* cm"-' 
and is accessible from the ground. Dense gas is also found in 
the circumstellar disk on scales of a few tens to hundreds AU, 
but such regions are generally diluted by an order of magnitude 
in single-dish observations. In contrast, molecular lines with 
much lower critical densities, such as the low excitation C"^0 
transitions, contain much higher contributions from low density 
material. This makes these lines well-suited as column density 
tracers for large-scale cloud material and the cold outer regions 
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Fig. 1. The L 1688 core in Oph iuchus. In grey-scale and contours the 850 //m SCUBA map as published by lJohnstone et al] (l2000h 
and iDi Francesco et all (l2008l) is shown. The locations of all the sources as observed in this study are shown, except for C2D- 
162527.6, Haro 1-4, C2D-162748.2 and IRS 63. C2D-16274.1 is marked as C2D-2 and C2D-162748.2 as C2D-3. The dashed line 
indicates the Oph ridge. 
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Table 1. Sample of embedded sources in L 1688, with additional known disks included. 



Source 


Other names 


Coordinates (J2000) 


Ref." 








RA 


Dec 




C2D-162527.6 


SSTc2d J162527.6-243648 


16:25:27.6 


-24:36:48.4 


2 


0.36 


GSS 26 




16:26:10.4 


-24:20:58 


1 


-0.46 


CRBR 2315.8-1700 




16:26:17.2 


-24:23:45.1 


2 


0.69 


CRBR 2317.3-1925.3 


SKS 1-10 


16:26:18.8 


-24:26:13 


1 


-0.56 


VSSG 1 


Elias 20 


16:26:18.9 


-24:28:22 


1 


-0.73 


GSS 30 


Elias 21/GY 6/GSS 30-IRSl 


16:26:21.4 


-24:23:04.1 


2 


1.46 


LFAM 1 


GSS 30-IRS3 


16:26:21.7 


-24 22 51.4 


2 


0.73 


CRBR 2324.1-1619 




16:26:25.5 


-24:23:01.6 


2 


0.87 


VLA 1623 


VLA 1623.4-2418 


16:26:26.4 


-24:24:30.3 


2 


no 


GY51 


VSSG 27 


16:26:30.5 


-24:22:59 


1 


0.05 


CRBR 2339.1-2032 


GY91 


16:26:40.5 


-24:27:14.3 


2 


0.45 


WL 12 


GY 111 


16:26:44.0 


-24:34:48 


1 


2.49 


WL2 


GY 128 


16:26:48.6 


-24:28:39 


1 


0.02 


LFAM 26 


CRBR 2403.7/GY 197 


16:27:05.3 


-24:36:29.8 




1.27 


WL 17 


GY 205 


16:27:07.0 


-24:38:16.0 


1 


0.61 


Elias 29 


WL15/GY214 


16:27:09.6 


-24:37:21.0 


1 


1.69 


GY 224 




16:27:11.4 


-24:40:46 


1 


-0.05 


WL 19 


GY 227 


16:27:11.9 


-24:38:31.0 


1 


-0.43 


WL20S 


GY240 


16:27:15.9 


-24:38:46 


1 


2.75 


IRS 37 


GY244 


16:27:17.6 


-24:28:58 


1 


0.25 


WL3 


GY249 


16:27:19.3 


-24:28:45 


1 


-0.03 


IRS 42 


GY 252 


16:27:21.6 


-24:41:42 


1 


-0.03 


WL6 


GY 254 


16:27:21.8 


-24:29:55 


1 


0.72 


GY 256 




16:27:22.0 


-24:29:39.9 




-0.05 


IRS 43 


GY265 


16:27:27.1 


-24:40:51 


1 


1.17 


IRS 44 


GY 269 


16:27:28.3 


-24:39:33.0 


1 


2.29 


Elias 32 


IRS 45/GY 273A'SSG 18 


16:27:28.6 


-24:27:19.8 




-0.03 


Elias 33 


IRS 47/GY 279A^SSG 17 


16:27:30.1 


-24:27:43 


1 


-0.12 


IRS 48 


GY 304 


16:27:37.2 


-24:30:34 


1 


0.88 


GY 312 




16:27:38.9 


-24:40:20.5 




0.64 


IRS 51 


GY3 15 


16:27:40.0 


-24:43:13 




-0.15 


C2D- 162741. 6 


SSTc2d 1162741.6-244645 


16:27:41.6 


-24:46:44.6 




0.32 


C2D- 162748.2 


SSTc2d 1162748.2-244225 


16:27:48.2 


-24:42:35.6 




1.55 


IRS 54 


GY 378 


16:27:51.7 


-24:31:46.0 




0.03 


IRAS 16285-2355 




16:28:21.6 


-24:36:23.7 


2 


1.23 


C2D- 162857.9 


SSTc2d 1162857.9-244055 


16:28:57.9 


-24:40:54.9 


2 


0.67 


IRS 63 




16:31:35.7 


-24:01:29.5 


2 


0.14 


Known Disks 


Haro 1-4 


DoAr 16 


16:25:10.5 


-23:19:14.5 


2 


-0.89 


DoAR 25 


GY 17 


16:26:24.0 


-24:43:09.0 


1 


-1.12 


OphE MM3 




16:27:05.9 


-24:37:08.2 


2 


-0.33 


SR21 


Elias 30 


16:27:10.2 


-24:19:16.0 


1 


-0.79 


CRBR 2422.8-3423.8 


CRBR 85 


16:27:24.8 


-24:41:03.0 


1 


1.01 


IRS 46 


GY274 


16:27:29.7 


-24:39:16.0 


1 


0.18 


SR9 


IRS 52/Elias 34 


16:27:40.5 


-24:22:07.0 


1 


-1.07 


2MASS 16282 




16:28:13.7 


-24:31:39.0 


2 


-1.55 



" Refer ences for po s itions : l lBontemps et alj ( l200ll) . 2: c2d delivery document 
* From lEvans et al.l ( l2009h . See also § 5.5 



of the protostellar envelope (e.g. |j0rgensen et al.l l2002h . Both 
the HCO^ and C"*0 data have velocity resolutions of 0.1 km 
s ' or better, thus allowing foreground clouds to be identified. 



Most sub-millimeter line data so far have been single 
pixel spectra toward the YSO with at best a few positions 
around specific YSOs. The recently commissioned HARP-B 
instrument is a 16-pixel receiver, operating in the 320 to 370 
GHz at mospheric window allowing rapid mapping of small (2') 
regions (ISmith et al.ll2003h . HARP-B is mounted on the James 



Clerk Maxwell Telescope (JCMTfl 

We present here HARP-B maps of all Class I sources in the 
L 1688 region in C'**0 3-2 and HCO+ 4-3. The combination of 
these two molecular lines allows us to differentiate between pro- 
tostellar envelopes, dense cores and foreground cloud material, 
as well as edge-on disks. The goal of this paper is to character- 
ize the envelopes of the embedded source population of L 1688, 
as well as present a new method for identifying truly embed- 



' The James Clerk Maxwell Telescope is operated by The Joint 
Astronomy Centre on behalf of the Science and Technology Facilities 
Council of the United Kingdom, the Netherlands Organisation for 
Scientific Research, and the National Research Council of Canada. 
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Fig. 2. SCUBA 850 fim continuum (left) and HCO+ 4-3 (right) 
integrated intensity maps plotted from bottom to top for GSS 30 
IRSl, VLA 1623, Elias 33 and IRS 43. The IR source positions 
are marked in the SCUBA maps with white crosses. Note that 
the absolute scale of these maps is a factor 4 higher than those 
in Fig. 3 and 4. The 850 /im continuum flux for VLA 1623 has 
been scaled down by a factor 4. In the map of GSS 30, LFAM 1 
is located north of the GSS 30 IRSl position. Elias 32 is located 
at the north-west side of the Elias 33 map. CRBR 2422.8-3423 is 
located to the south-west of IRS 43. All four main targets show 
strong HCO^ and SCUBA peaking within 5" of the IR position 
at (0,0), characteristic of truly embedded sources. 



ded sources and separate them from (obscured) edge-on disks 
using the dense gas present in embedded YSOs. In § 2, a sam- 
ple of Class I sources in the L 1688 core is selected. In § 3, we 
discuss the details of the heterodyne observations carried out at 



flie JCMT and the Atacama Pathfinder Experiment (APEXfl as 
well as the supplementary observations obtained in a continuum 

1. e. wideband mode. § 4 presents the maps and spectra and in § 
5 we analyze the properties of the gas and dust of the sources in 
the sample. The environment around the YSOs, column density, 
envelope gas and concentration of the HCO^ are discussed. In 
§ 6, we present a new method for identifying embedded YSOs 
from (edge-on) disks and apply this method to the sample. This 
classification is then compared to traditional methods as well as 
other recently proposed methods. The main conclusions of the 
paper are given in § 7. 

2. Sample Selection 

Of the known (embedded) YSO population within L 1688, 45 
objects were selected for our sample using several criteria. First, 
we require all sources be located within the 850 jum dust con- 
tinuum map made by the COMPLETE project of L 1688 using 
SCUBA on the JCMT ( John stone et al. 2000). 

Second, sources must be included in the area covered 
by the c2d program using IRAC and MIPS on Spitzer 
dEvans etalj r2003. Padgett etal. 20 0i). All so u rces classified 
as Class I in either Andre & Mo ntmerld d 19941) . iBarsonv et al.l 



( 19971). iBontemps et al. (2001) or the c2d delivery document 



(E vans etal J l2007lPI are included with luminosities > 0.04 Lq. 
Although these 41 objects have been classified as Class I in one 
or more of these papers, only 6 of these have been classified 
consistently as Class I in all studies. Most other sources are clas- 
sified as either Class II or Flat spectrum sources at least once. All 
such sources are included in the analysis of this paper, but con- 
versely, our sample does not include all Flat spectrum or Class 
II sources listed in the c2d survey. Known edge-on disks such as 
2MASS 16282, IRS 46, OphE MM3 and CRBR 2422.2-3423, 
are among these 41 sources and are included within the sample 
to illustrate the results of our method for such sources. The sam- 
ple should not contain any reddened background main-sequence 
stars which are readily identified in the c2d analysis. However, 
other background sources, such as AGB stars or background in- 
frared galaxies, may be present. 

Four sources with a2-24f/m> 0.0 were found only in the 
recent Spitzer observations. These are SSTc2d J162527.6- 
243648, SSTc2dJ162741.6-244645,SSTc2dJ162748.2-244225 
and SSTc2d J162857.9-24405. The names C2D-162527.6, C2D- 
162741.6 , C2D-162748.2 and C2D-162857.9 are adopted. VLA 
1623-2418 is included as an embedded Class 0, but it is gen- 
erally absent from the above studies. In addition, IRS 63 was 
included, although it is not part of the L 1688 core, due to its 
interesting characteristics as a Class I source as observed with 
the SubMillimeter Array (SMA). These recent interferometric 
results suggest that this source is an embedded YSO with lit- 
tle envelope material left, and, as such, presents an interesting 
test case for the proposed met hod of identifying truly embedded 
sources dLommen et al.ll2008l) . Four disk sources within or near 
to the L 1688 region, Haro 1-4, DoAr 25, SR 9 and SR 21, were 
included to serve as a sample of known disk sources. 

The final source sample, which covers all potential Class 
I s ources, can be found in Table [1] The positions as found 
in iBontemps et alj d200Ih were used where available. If the 



- This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fur Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 



' ,http://ssc. spitzer.caltech.edu/, 



6 T.A. van Kempen et al.: The Class I population in the Ophiuchus clouds 



SCUBA C'*tl3-2 




Jyiitwni K km f " 4C tern 

^^^^^H^TJ ^^H^I^TH ^^^^^^KZD 

t Uf BED it 4 I 



Fig. 3. SCUBA 850 fj.m continuum (left), HCO+ 4-3 (middle) and C'^O 3-2 (ng/if) integrated intensity maps plotted for (from 
bottom to top) IRS 44, LFAM 26, WL 12, IRAS 16285, IRS 54 and IRS 37. The locations of the IR sources ai-e indicated in the 
SCUBA maps with white crosses. In the IRS 44 map, IRS 44 is indicated with a black cross, and IRS 46 with a white one. In the 
LFAM 26 field, OphE MM3 is located south of LFAM 26. In the IRS 37 field, WL 3 is located west of IRS 37. The sources at (0,0) 
in this figure show weaker HCO^ images than in Fig.|2] but are still peaking on the IR positions in both SCUBA and HCO^, and 
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Fig. 4. SCUBA 850 ^im continuum (left), HCO+ 4-3 (middle) and C'^^O 3-2 (Wg/if) plotted for (from bottom to top) IRS 48, IRS 
63, IRS 42, WL 19, IRS 51 and WL17. The locations of the IR sources are indicated in the SCUBA maps with white crosses. The 
embedded source Elias 29 is shown with a black cross in the field of WL 17. The peak 850 fim flux density for IRS 63 is 1.1 Jy 
beam"'; it was not scaled to accentuate the extended dust emission. The sources at (0,0) in this image are mostly unresolved in 
SCUBA and show little or no HCO^ emission related to the source. They are either transitional or confused sources or (edge-on) 
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fusion exists due to nearby I R sources, the pos ition as found 
in the c2d delivery document ("Evans et al."2007') was adopted. 
For sources in common, most positions agree within 3". Fig. [T] 
shows the distribution of the sample as plotted on the SCUBA 
850 ^m map of L 1688 dJohnstone et aljioOOl lRidge et al.l2006l 
iDi Francesco et al.ll2008l) ! 
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Fig. 5. SCUBA 850 //m continuum (left) and HCO+ 4-3 (ng/if) 
plotted for (from bottom to top) CRBR 2339.1, WL 2, C2D 
162527.6, GY 312, C2D 162741.6 and CRBR 2315.8. The lo- 
cations of the sources are indicated in the SCUBA maps with 
white crosses. These sources have either unresolved or little or 
no emission in the SCUBA 850 yum and HCO+ 4-3 maps at (0,0), 
making them disk sources in our classification. 



source was not included in iBontemps et al.l (1200 lb . or if con- 



Fig. 6. HCO^ 4-3 spectra extracted from the HARP-B maps in 
a 15" around the IR positions. Sources classified as embedded 
in § 6.2 are shown in the left column. In the right column, the 
sources classified as confused or disks are plotted. 

The table includes the spectral index a2-24fim given in 
lEvans et al. (2009), calculated from the 2MASS, IRAC and 
MIPS (24 jum only) fluxes. No qtir could be determined from 
Spitzer data for VLA 1623 (too faint) 
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3. Observations 

3.1. Gas line maps 

The majority of the sample was observed in HCO^ 4—3 
(356.7341 GHz) and C"*0 3-2 (329.3305 GHz) using the re- 
cently commissioned 16-pixel heterodyne array receiver HARP- 
B on the James Clerk Maxwell Telescope (JCMT). The high 
spectral resolution mode of 0.05 km s ' available with the 
ACSIS back-end was used to disentangle foreground material 
as well as any contributions from outflowing material. Spectra 
were subsequently binned to 0. 15 km s"' . 

HARP-B observations of 30 sources in 21 fields were carried 
out during July and August 2007 in C'**0 3-2 and HCO+ 4-3 
under weather conditions with an atmospheric optical depth, 
T225GHZ ranging from 0.035 to 0.08 (precipitable water vapor of 
0.7 to 1.6 mm). The fields were observed down to a rms noise 
of 0.1 K in a 0.5 km s ' bin. The HARP-B pixels have typical 
single side-band system temperatures of 300-350 K. The 16 
receivers are arranged in a 4x4 pattern, separated by 30". This 
gives a total foot print of 2' with a spatial resolution of 15", the 
beam of the JCMT at 345 GHz. The 2' x 2' fields were mapped 
using the specifically designed jiggle mode HARP^. A position 
switch of typically 30' in azimuth was used, with larger throws 
if needed. The Class source IRAS 16293-2422 was used as 
a line calibrator and pointing source. Calibration errors are 
expected to dominate the flux uncertainties, and are estimated at 
20%. Pointing was checked every two hours and was generally 
found to be within 2-3". The map was re-sampled with a pixel 
size of 5", which is significantly larger than the pointing error. 
The main-beam efficiency was taken to be 0.67. Data were 
reduced using the STAI^INK package GAIA and the CLASS 
reduction package. 



3.2. Gas single pixel spectra 

Supplementary data of HCO^ 4-3 were taken at the APEX tele- 
scope during July 2007, using the APEX-2a receiver All sources 
were observed for which no HARP-B data were taken, except 
GSS 26 and GY 5 1 . The APEX observations were done in excel- 
lent weather conditions with T225ghz ranging from 0.01 to 0.04 
(PWV 0.2 to 0.8 mm). Single spectra were taken with a spectral 
resolution of 0.4 km s"' down to an rms of 0.3 K. Calibration er- 
rors are estimated to be ~ 20%. The APEX 12-m dish is slightly 
smaller than the JCMT 15-m dish, producing a beam of 18" in- 
stead of 15". Pointing errors are <4". Beam efficiency is 0.70. 

For the sources GSS 30, VLA 1623, WL 20S, IRS 42, IRS 
43, Elias 32 and Elias 33, C'^^O 3-2 spectra were obtained from 
the CADC archiv^ These spectra were taken with the RxB re- 
ceiver during September 2005. 

3.3. Dust maps 

The 850 jum continuum data of the Ophiuchus region, ob- 
tained within the scope of the COMPLETE project using the 
SCUBA instrument on the JCMT, were used to characterize 
the dust in the environments around all sources of the sample 
dJohnstone et all l2000i iRidge et all l2006i iDi Francesco et all 
120081). The most recent map, version 3, was used to extract 
the information (see Fig. [1). This version includes a correction 
for the chopped out-emission. Integrated fluxes within regions 



with radii ranging from 25" to 40" were extracted from the 
map (see Table |9]l. The exact radii were selected by calculating 
the FWHM to the peak flux of e ach source. Com parison with 
the map and published fluxes by IJohnstone et al.l (|2000|) within 
similar radii found by a clumpfind routine showed that fluxes 
agreed within 5%, significantly better than the calibration error 
of 20% for SCUBA. The sensitivity of the SCUBA map is such 
that sources down to 90 mJy (3cr) can be detected. 

In addition, 11 sources were observed using the SHARC- 
II 350 i-im continuum instrument at the Caltech Submillimeter 
Observatory (CSOf|during April 2003 (iDoweU et al.l2003b . The 
array has 12x32 pixels, spanning 4.85" per pixel, which results 
in a footprint of 1' x 2.6' and a beam size of 9". The sources ob- 
served with the CSO were GSS 30, WL 12, Elias 29, VLA 1623, 
OphE MM3, WL 20S, WL 6, IRS 43, VSSG 17, IRS 5 1, 2MASS 
16282. The data were reduced with the CRUSH package and cal- 
ibrated using observations of Mars and Saturn. If no planets were 
available during the specific night, IRAS 16293-2422 was used 
to calibrate the data. The calibration errors are estimated to be 
on the order of 30^0% and assumed to dominate the error in 
the flux estimate over other instrumental errors or S /N consider- 
ations. The sensivity of the SHARC-II maps is such that sources 
down to 400 mJy (3cr) can be detected. For a typical dust flux 
law oc v^*^, this implies that the SHARC-II data are a factor of 4 
more sensitive to low-mass sources than the SCUBA data. 



3.4. SED and IRS spectra 

The L 1688 core has been targeted by a large number of 
continuum surveys, covering wavelengths from 2 to 1300 
lum. 2MASS (1.25, 1.66 and 2.2 fim) and Spitzer-IRAC (3.6, 

4.5, 5.8 and 8 jum) and Spitzer-MIPS (2 4 ^m) fluxes wer e 
obtained from the c2d delive ry document jEvans et al.ll2007l) . 
Recendy, iPadgett et all (|2008|) presented Spitzer-MIPS 70 and 
160 fim fluxes. The MIPS 160 jum data were not used in our 
analysis since a large part of the map was confused by either 
striping or saturation. For the MIPS 70 jum data, we chose 
to go back to the original data and retrieve the 70 fim fluxes 
source-by-source in a typical 10" radius using PSF photometry 
without the extended contribution subtracted. This avoids 
potential errors in the flux estimates from automatic extraction 
in crowded or confused regions such as the Ophiuchus core. In 
addition to the SCUBA (850 yum) and/or SHARC-II (350 fim) 
fluxes, a limited number of sources was also observed using 
SHARC-II and/or SCUBA at 450 fim in a recent survey of disks 
in the L 1688 core, with flu xes given in the central 15" beam 
(lAndrews & Williamsl l2007h . Where available, archival data 
at millimeter wavelengths were used from SEST and IRAM 
30m, as rep orted i n Andrews & Wi lliam s (2007) and originally 
obtained by lAndre & Montmerl d (Il994 . Jensen et al. ( 199^, 
iNue mberger et al.l (il99&) . .Motte et al.l ([T998i) and iStanke et al] 
(200^^ 



A number of sources in our sample have also been observed 
with various observing programs using the IRS instrument on 
Spitzer at 5-40 //m (PID= 2, 172 and 179). These spectra are 
included in this study to confirm the continuum fluxes in the 
wavelength range of IRS, and to investigate the presence of ice 
absorption and silicate absorption or emission. 



See JCMT website http://www.jcmt.jach.hawaii.edu/ 
' See ,http://www.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/jcmt/1 



http://www.submm.caltech.edu/cso/ 
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Fig. 7. C'^O 3-2 spectra extracted from the HARP-B maps at the source positions. Spectra are arranged in random order to illustrate 
the independence of C"^0 to the environment instead of the protostellar envelope. 
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Fig. 8. HCO^ 4-3 spectra at the source positions obtained at APEX. Elias 29 and IRS 43, with clear detections, are the only sources 
found to be embedded in § 6.2. The detection seen for CRBR 2317 is attributed to the cloud material and is classified as a disk. 



4. Results 

4.1. Gas maps 

Fig. |2] to ID present the 21 observed HARP-B fields for HCO+ 
4-3 and C'-^O 3-2, together with the SCUBA 850 /im continuum 
map. Note that the sources in Fig. |2] are much brighter and 
have difi'erent intensity scales. The integrated intensity and peak 
temperatures for all sources are listed in Table |2l columns 2 to 

5. Column 7 states in which figure(s) the sources are mapped 
or in which figure(s) its spectra are plotted. Fig. |6](HCO^) and 
|7](C'^0) show spectra extracted from the maps at the positions 
of the sources. About 70% of the sample was observed in C'^O 
3-2. The sources for which APEX-2a HCO^ spectra were taken 
can be found in Fig.[8l 



Using the molecular line emission maps, the sample can be 
divided roughly into three groups. In the first (see Fig. |2] and 
[3]l, HCO^ peaks at or close to the source positions, with peak 
intenstities up to a few K. For these sources, GSS 30 IRSl, 
LFAM 1, VLA 1623, WL 12, LFAM 26, IRS 37, WL 3, WL 
17, EHas 33, Elias 32, IRS 44, IRAS 16285 and IRS 54, HCO+ 
seems to be extended on scales of a few 10". C'^O is always 
present in these maps, and also extended on scales of a few 10" 
to 1'. 

Note that the HCO^ maps in some case show minor off- 
sets w.r.t. the IR positions. For offsets < 8" this can be at- 
tributed to either the pointing accurac y of the JCMT or the 
differences between IR positions from IBontemps et alj (l200lh 
and c2d. For a few sources, in particular VLA 1623, IRS 44, 
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Fig. 9. Dust maps at 850 yLun extracted from the JCMT-SCUBA COMPLETE map of GSS 26, Ellas 29 and IRS 43; CRBR 2317.3- 
1925, VSSG 1 and GY 224; WL 6, C2D-162748.2 (C2D-3) and C2D-162857.9 (C2D-4); GY 51, SR 21 and SR 9; Haro 1-4, DoAr 
25 and 2MASS 16282. Note that the emission toward GY 51 is scaled down by a factor of 4. 
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Table 2. Results from HARP-B (HCO+ 4-3 and C'^^O 3-2), APEX-2a (HCO+ 4-3) and RxB (C^^O 3-2) observations at the source 
position. 



Source 


HCO+ 4-3" 




C"*0 3-2 




Notes 




/ Tmb dV 


Tmb 


/ Tub dV 


Tmb 






[Kkms-'] 


[K] 


[Kkms-'] 


[K] 




C2D-162527.6 


- 


<0.09 


- 


- 


SeeFig.|5| 


GSS 26 


- 


- 


- 


- 


See Fig. [9] 


CRBR 2315.8-1700 


6.3 


5.2* 


- 


- 


See Fig|5] 


CRBR 2317.3-1925 


1.5* 


2.0'' 


- 


- 


See Fig. [9] 


VSSG 1 


- 


<0.4 


- 


- 


See Fig. [9] 


GSS 30 


16.7 


9.1 


10.7'^ 


6.3" 


See Figl2]and[i0] 


LFAM 1 


18.4 


10.8 


- 


- 


See Fig|2]and[i0l GSS 30 maps 


CRBR 2324.1-1619 


- 


- 


14.5 


12.2 


See Figl2]and[i0l VLA1623 maps 


VLA 1623 


16.2 


9.0 


10.4 


17.9'' 


See Figl2]and[T0] 


GY51 


- 


- 


- 


- 


See Fig. [9] 


CRBR 2339.1-2032 


0.55 


0.45* 


- 


- 


See FiglU 


WL 12 


0.87 


1.4 


4.8 


3.9 


See Figl3]and[iO] 


WL2 


0.3 


0.18* 


- 


- 


See FigIS] 


LFAM 26 


1.95 


1.8 


6.4 


3.3 


See Fig|3] 


WL 17 


0.6 


0.45 


- 


- 


See Figil 


Elias 29 


4.5'-" 


2* 


10.4" 


4.3" 


See FigBl WL17 map and Fig.|9]and[l0] 


GY 224 


- 


<0.27'' 


- 


- 


See Fig. [9] 


WL19 


- 


<0.1 


6.3 


4.5 


See Figil 


WL20S 


- 


<0.26'' 


5.0" 


2.3" 


See Fig. [H 


IRS 37 


3.7 


2.5 


9.55 


5.7 


See FiglU 


WL3 


3.7 


2.2 


8.8 


6.3 


See FigUl IRS 37 map 


IRS 42 


0.45 


0.7* 


4.0 


1.2* 


See Figia 


WL6 


0.92 


0.6 


5.25 


4.5* 


SeeFig.[9]and[T0] 


GY 256 


0.75 


0.45* 


5.4 


4.5* 


See Fig.[9]and[T0l WL6 maps 


IRS 43 


4.5'' 


1.6" 


6.5" 


3.2" 


SeeFig.[2l9land[70l 


IRS 44 


3.4'' 


0.5'' 


5.8 


2.8 


See Fig|3] 


Elias 32 


5.6 


4.3 


3.1" 


4.5" 


See Fig|2]and[T0l Elias 33 maps 


Elias 33 


6.3 


5.7 


6.0" 


3.0" 


See Figl2]and[i0] 


IRS 48 




<0.09 


2.6 


2.3 


See Figil 


GY312 




<0.1 






See FigIS] 


IRS 51 


0.75 


0.75 


2.0 


2.8* 


See Figi]and[TO] 


C2D-162741.6 




<0.1 






See FigIS] 


C2D-162748.1 




< 0.34'' 






See Fig. [9] 


IRS 54 


0.53 


0.45 


3.4 


3.3 


See FigE] 


IRAS 16285-2355 


2.1 


3.0 


2.3 


3.0 


See FigH] 


C2D-162857.9 




< 0.34'' 






See Fig. [9] 


IRS 63 


0.75 


1.2 


1.7 


3.3 


See Figil 


Disks 


Haro 1-4 




<0.29'' 






SeeFig.|9| 


DoAR 25 




<0.28'' 






See Fig. [9] 


OphE MM3 


1.7 


1.6* 


5.5 


3.0* 


See FigUl LFAM 26 map 


SR21 




<0.28'' 






See Fig. [9] 


CRBR 2422.8-3423 


1.9 


1.4* 






See Fig. [2] IRS 43 map, and FigilNW comer of IRS 42 


IRS 46 




<0.09 


2.9 


2.1 


See FigUl IRS 44 map 


SR9 




<0.24'' 






See Fig. [9] 


2Mass 16282 




<0.35'' 






SeeFig.[9]and[T0] 



" Intensities marked with a * do not peak at source position 

APEX-2a receiver. Upper limits (2o") in 0.4 km s"' bin 
" JCMT RxB data 

'' Outflowing gas detected (width ~20 km s"') 

" C'^O data taken from .Boogert et al.. C20d3) , obtained with the CSO. 



larger offsets up to 12" were found. These can be accounted 
for by potentially three effects. First, self-absorption of HCO^. 
Second, these sources show strong outflow emission, which 
peaks off-source and significantly influences the spectrally in- 
tegrated maps. Third, molecular emission can be suppressed by 
the freeze-out of molecules in the cold part of the cloud. 



A second group (Fig. |4]and|5]l lacks detections in HCO^, 
down to a Hmit 0.1 K. The sources, C2D-162527.6, WL 19, 
IRS 46, IRS 48, GY 312 and C2D-162741.6 belong to this 
group. The third group (Fig.|4]and|5]i of sources contains CRBR 
2315.8-1700, CRBR 2339.1-2032, WL 2, OphE MM3, IRS 
51 and IRS 42. These show HCO^ detections at the source 
position. The HCO^ is extended, but there is no sign of a peak at 
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Fig. 10. Dust maps at 350 //m observed with the SHARC-II instrument at the CSO. In white, the contours for the JCMT-SCUBA 
850 jjm data are over-plotted, normalized for the maximum emission in the map in contours of 10%, 30%, 50% and 70 %. The 
shape of the contour levels at 850 yum agrees well with those found at 350jum. 



the source positions. For IRS 63 and WL 6, detections of HCO^ 
are marginal at ~3 cr. 

The HCO^ 4-3 integrated intensities J TMBdV in protostel- 
lar envelopes range from as high as 18.4 K km s ' for LFAM 1 
to 0.75 K km s ^ for IRS 63. In § 5.3 the HCO+ 4-3 integrated 
intensities are compared with the column densities derived from 
the C'^O 3-2. This shows that most embedded sources with 
strong HCO^ emission are located in regions with high column 



densities. The sole exception is IRAS 16285-2355, which is not 
located in a region with high column density (see Fig. 1 andO. 
However, it is possible to find sources with little or no HCO^ 
4-3 emission, such as WL 19, in regions with similarly high 
column density. 
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4.2. Dust maps 

For sources for which molecular emission was taken with the 
HARP-B array, the spatial extent of the dust is shown in the left 
column in Fig. |2l|5] The sources for which no HARP-B maps 
were obtained can be found in Fig. |9] The SHARC-11 maps are 
shown in Fig. [TO] In these maps, the white contours show the 
850 SCUBA emission. The integrated fluxes in a 25" beam are 
given in Tables [3] and |9l 



Table 3. Source fluxes at 350 fim in a 25" beam observed with 
SHARC-11 at the CSO. 



Source 


Flux 


Concentration" 




(Jy) 




GSS 30 


2.5 


0.58 


WL12 


6.7 


0.59 


Ellas 29 


4.3 


0.72 


VLA1623 


40 


0.8 


WL20S 


0.82 


0.64 


WL6 


<0.4 


U 


IRS 43 


2.1 


0.62 


Elias 33 


4.4 


0.55 


IRS 51 


1.4 


0.64 


2MASS 16282 


<0.2 


U 



" See §5.1; U=not determined since only flux upper limits obtained. 



Both the SHARC-II and SCUBA images clearly show that 
the dust is extended on scales of at least a few arcminutes at 
many positions in Ophiuchus. The smaller beam and higher 
frequency of the SHARC-11 observations is able to both resolve 
smaller envelopes (e.g., Elias 33) or confirm that other sources 
(e.g., WL 20S) are not resolved at 350 jjm down to 9". The 
extended emission originates from co l d dust in the parenta l 
clo ud, as mapped bylMotte et aP (Il998h . iJohnstone etail (l2000h 
and IStanke et al.l ((2006'). This cloud material can exist close to 
or in the line of sight of many of the Class 1 sources, but is not 
necessarily associated with a protostellar envelope, as is often 
assumed. A good example is the source WL 19 (Fig.|4|i, where 
dust emission is found close to the position of the source, but 
does not peak at the position of WL 19. This emission comes 
from a prestellar core within the Ophiuchus ridge and not a 
protostellar envelope. 

Based on the spatial extent of the dust, the sample can be 
divided into four groups. First there are sources with spatially 
extended dust emission profiles, peaking at the source position. 
Examples ai-e GSS 30 IRSl, Elias 29 and IRS 43. A second 
group shows extended dust emission, but with no central peak 
at the position of the source. A good example is GY 51 (see Fig. 
|9ll, located close to the Oph-A core. Some Class 1 sources only 
show unresolved dust emission at the source position, such as 
IRS 51 (see FigfTOll. A final group shows no emission at 350 or 
850 fim down to our sensitivity limits. Good examples is Haro 
1-4. 

Cloud material can contribute significantly or even dominate 
the emission on the scales probed by SCUBA and SHARC-11. 
Two fluxes at 850 /vm are therefore given in Table |9] The first 
flux is within the 15" beam obtained from Andrews & Williams 
(l2007h . If no value was given by them, it was extracted 
from the COMPLETE map for the central beam. Our ex- 



tracted fluxes agree t o with in 15% for sources also listed in 
[Andrews & Williamsl (l2007l) . even though we include estimates 
to reconstruct the chopped-out extended emission. The second 
flux is extracted from the COMPLETE map for larger apertures 
encompassing the envelopes, up to a 25" (~3,000 AU) radius, 
the typical envelope extent where the temperature and density 
drop to that of the surrounding cloud. The first number should 
thus be considered a lower limit on the total envelope and disk 
emission, while the second is an equivalent upper limit. For 
disk sources, all emission originating from the source is located 
within 15" (<1,800AU). 



4.3. SED and IRS spectra 

Table |9] summarizes all the Spitzer and (sub-)millimeter fluxes. 
The 2 fj.m to 1.3 mm flux densities for all 45 sources are plot- 
ted in Fig. [TT] In addition, low-resolution IRS spectra are over- 
plotted between ~5-30 fim, where available. In the on-line ap- 
pendix, blow-ups of the IRS spectra are presented. 

5. Analysis 

5.1. Concentration 

One method to distinguish embedded YSOs from edge-on disks, 
starless cores or background sources such as AGB stars and 
galaxies, is to use a measure of how centrally concentrated the 
emission is. This meth od has been used previously in the anal- 
ysis of SCUBA da ta (John stone etal.ll200lL IWalawender et al.1 
l2005l |j0rgensen et al...2007.) by defining a concentration param- 
eter: 

c = i-Hl^ (1) 

where B is the beam-size, ^850 is the total flux at 850 yum 
within the envelope, R^^^ is the observed radius of the envelope 
and F() is the peak flux of the envelope. The observed radii 
were selected by calculating the FWHM to the peak flux of 
each source. The level of concentration potentially provides an 
indication of the evolutionary stage. Highly concentrated cores 
(C > 0.75) have a very high probability to contain protostars, 
while cores with low concentration (C < 0.4) are often starless, 
as indicated by the absence of near-IR emiss ion. The observed 
radii in the 850 fim maps were defined by iKirk et aT\ (l2006h 
as = (A/tt)'''^, with A the area of the core within the lowest 
contour of their clumpfind routine. The SHARC-II maps also 
allow a concentration at 350 ;um to be made, as seen in Table [3] 
These concentrations are found to be within 10% to those found 
for the 850 //m data, with the exception of WL12 which is more 
than a factor 2 higher in the SHARC-II 350 yum data. 

A similar concentration parameter can also be calculated for 
the distribution of HCO^, using integrated intensities between 
Vlsr - 5 to +5 km s"' : 

Chco. = 1 - (2) 

with 5hco+ the spatially and spectrally integrated intensity 
/ TundV of HCO+ over the entire envelope with radius Robs, 
and Tmb the peak intensity in the central beam. Since no 
clumpfind routine was used to characterize the spectral maps, 
the observed radii of the HCO^ concentration were set equal to 



T.A. van Kempen et al.: The Class 1 population in the Ophiuchus clouds 



15 



CN 
I 

E 

o 



CD 
Ll 



1 * 


LFAM 1 

1 


Vl Al fi9 

X 

X 

* X 

*. 1 


WL 1 2 

X 

X 

X 
X 

X 

X 

1 * 


1 FAM9fi 
i_ r /A IVI t u 

X 1 * 


WL 1 7 

X 

1(T) 


^ Fl IA*^9Q 

* * 

X 

X 

1 5' 


IRS37 

5K 

X 

1 


WL3 

1(T) 


WL6 

X 

X 
X 

* 1 X 

X 


IRS43 

X 

X « 
X 1 * 


IRS44 

' X 

« X 
X 

1 X 


Fl IA*^^9 

X 

x^^ 

X X 
X 

1 


Fl lA*^^^ 

X^^ ^ 
X 

X * 

X 

1 


IRS54 


IRA^^I fi9RS 

X 

" 1 


IRS63 

X 

X * 

X 

X 

X 

HT) 


rRRR9^94 1 

X 

c 


GY51 

JK X 
X 

X 

X 

c 

X 


WL 1 9 

X * 


IRS42 

X 

X 

X 

X 

* X 

c 


GY256 


r9ni fi9S97 fi 

X 

2 


GSS26 
2 


rRRR9'^1 S R 

* 2 * 


rRRR9'^1 7 

X 

2 


VSSGl 

X 

X 

2 


rRRR9'^'^Q 1 

X X, 

2 


Wl 9 


X 

X 

* 9 


Wl 90*^ 

./^'^ 

3IW X 

X 

X 

" 2 


IRS48 

^J^V' X 
X 

X 

2 X 


1 9 

* X 

X 
X 


IRS51 

X 

X 

X 

X 

2 * 


r9nifi974i fi 

i(i«X X * 

X 

.X X 

2 


r9ni fi974R 9 

« 2 * 


X « 
X 

X 2 


2 X 


nnAR9S 
SI * 

*x 

2 


~ n c. ivi IVI J 

2 

X 

X 


'^R91 

>^%«/" 

X 

X 

2 * 


rRRR94.99 R 
f\ D r\ z t z z ,o 

X 

X 


IRS46 

2 


SR9 

2 

X 


9M1 fi9ft9 

£. IVI 1 \J£.KJ£. 

2 

X * 

« X 





1 10 10^ 10^ 10 10^ 10^ 10 10^ 10^ 



A [/xm] 



Fig. 11. The spectral energy distributions of all sources, supplemented by spectra taken with IRS. Our classification can be found 
near the center and bottom of each SED. 1 = Embedded Stage 1, 1(T) - Transitional from embedded to pre-main sequence, C = 
confused and 2 - (edge-on) disk. 



the FWHM of each spectrally integrated core map. Although 
different methods to define radii result in different values, tests 
showed that the corresponding changes in concentration factor 
are less than 5 %. 

Tables[3]and|4]list the concentrations of both the 850 jim data 
SCUBA maps and the HCO+ HARP-B maps, where available. 
Although the number of sources with a concentration parameter 
determined by SCUBA maps is greater than those for HCO^, 



it is clear that HCO^ is more concentrated than the continuum. 
Fig. [12] shows the distribution of the concentration parameters 
for spatially resolved sources. 

5.2. Environment 

In contrast with deeply embedded Class YSOs, the environ- 
ment around the embedded Class I YSOs has a large influence 
on the analysis of the generally weaker source data. Fig. 7 in 
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Table 4. Source properties derived from the gas and dust. The table has been rearranged compared with Table [1] to reflect the 
classification in §6.2 



Source 




N" 

(10"22 cm-2) 




Mass"-' 
(10-2 Mo) 


(10-' Lo) 


(K) 


r d 

<-850 


<-'HCO+ 


Embedded sources 


V.J O O ~J\J 


1.46 


23.3/19.4* 


2.6 


20.5 


33 


123 


0.55 


0.69 


LFAM r 


0.73 


50* 


2.6 


17.1 


8.3 


86 


0.59 


0.69 


VI A 1673 




441* 




111 


2.7 


12 


0.7 


0.72 


WL 12 


2.49 


10 5/15 2* 


3.1 


4.6 


34 


155 


0.26 


0.81 


LFAM 26 


1.27 


14.0/18.5* 


4.4 


4.5 


0.44 


238 


0.40 


0.69 


WL 17 


0.61 


8.4* 


1.0 


4.0 


6.7 


323 


0.22 


0.71 


Elias 29 


1.69 


22.7/33.6* 


1.5 


6.2 


25 


424 


0.49 


./ 


IRS 37'- 


0.25 


20.8 


1.0 


1.2 


3.8 


243 


0.58 


0.72 


WL3'' 


-0.03 


19.2/9.7* 


2.3 


2.9 


4.6 


192 


0.25 


0.67 


WL6 


0.72 


11.4 




<0.4 


8.5 


394 


U 


0.77 


IRS 43 


1.17 


14.2/20.7* 


2.3 


17.1 


10 


134 


0.58 


0.79 


IRS 44 


2.29 


12.6/30.2D* 


2.9 


8.0 


11 


140 


0.33 


0.80 


Elias 32 


-0.03 


6.8/140* 


3.7 


41.9 


5 


321 


0.37 


0.70 


Elias 33 


-0.12 


13.1/96* 


3.1 


28.2 


12 


460 


0.6 


0.64 


IRS 54 


0.03 


7.4 


1.2 


3.1 


7.8 


486 


0.8 


0.80 


IRAS 16285-2355 


1.23 


5.0/41.9D* 


3.1 


10.3 


0.57 


77 


0.45 


0.80 


IRS 63 


0.14 


3.7/84.0D* 




16.4 


13 


363 


0.75 


0.60 


Confused sources 


CRBR 2324.1-1619 


0.87 


31.6/70.5* 


5.2 


39.3 


0.44 


283 


U 


U 


GY51 


0.05 


27.7* 


4.5 


11.1 


2.3 


706 


U 


u 


WL19 


-0.43 


13.7/16.8* 


4.9 


5.5 


1.5 


730 


0.15 


0.53 


IRS 42 


-0.03 


8.7/14.2* 


1.3 


4.7 


14 


540 


0.44 


0.60 


GY 256 


-0.05 


11.8 




<0.2 


1.1 


104 


U 


0.70 


Disks 


C2D-162527.6 


0.36 






<0.4 


0.16 


1051 


u 


U 


GSS 26 


-0.46 


24.3D* 


1.8 


10.2 


3.2 


877 


0.65 


u 


CRBR 2315.8-1700 


0.69 


17.6* 


1.0 


5.1 


0.93 


456 


0.38 


0.74 


CRBR 2317.3-1925 


-0.56 


12.6* 


0.3 


1.1 


0.11 


1105 


0.85 


S 


VSSG 1 


-0.73 


20.1* 




2.0 


5.8 


930 


0.74 


SU 


CRBR 2339.1-2032 


0.45 


6.7* 




2.0 


0.55 


426 


0.55 


0.58 


WL2 


0.02 






1.0 


0.88 


573 


0.3 


0.58 


GY 224 


-0.05 


10.9* 


1.0 


2.4 


2.1 


614 


0.63 


SU 


WL20S 


2.75 


10.9/15.3D* 


1.0 


3.9 


4.5 


85 


0.57 


SU 


IRS 48 


0.88 


5.7/15.2D* 


2.1 


1.2 


81 


238 


0.70 


u 


GY312 


0.64 






<0.2 


0.73 


285 


U 


0.3 


IRS 51 


-0.15 


4.4/1 8.4D* 


2.3 


5.0 


8 


548 


0.27 


0.59 


C2D-162741.6 


0.32 






2.3 


0.28 


279 


U 


U 


C2D-162748.2 


1.55 






0.4 


0.17 


241 


U 


SU 


C2D-162857.9 


0.67 






5.6 


0.22 


402 


u 


SU 


Haro 1-4 


-0.89 






0.45 


5.6 


762 


u 


SU 


DoAR 25 


-1.12 


38.6* 




3.4 


4.7 


1422 


0.90 


SU 


OphE MM3 


-0.33 


12.0/13.4* 


3.4 


2.3 


0.17 


45 


u 


0.68 


SR21 


-0.79 


33.6* 




3.1 


1.8 


1070 


0.89 


SU 


CRBR 2422.8-3423 


1.01 


106.6* 


2.3 


20.5 


4.4 


157 


U 


./ 


IRS 46 


0.18 


6.3/15.1* 


2.8 


3.3 


1.9 


684 


U 


SU 


SR9 


-1.07 






<0.4 


8.7 


1588 


U 


SU 


2MASS 16282 


-1.55 






1.7 


0.84 


66 


U 


SU 



" Column densities marked with * are calculated from the SCUBA dust maps. Embedded sources where a large disk contribution is likely are 
indicated with a D. Note that no column densities could be calculated from the dust below 1x10"'' cm-^ . 

' Aa is an upper limit to change in 2-24 pm spectral slope due to the reddening by foreground clouds at this source. This is derived by considering 
the mean column density at 45" off the source position. 

Upper limit derived from the 850 pm data using the relation from lShirlev et al.l i2000h . with the assumption that all 850 pm emission originates 
within a protostellar envelope. 

U=Undetected or Unobserved, S=single position spectrum available only 

Circumbinary envelope. 

^ Although both Elias 29 and CRBR 2422.8-3423 were observed, these sources are located near the edge of a field and do not have sufficient 
coverage to obtain a concentration parameter for the HCO^. 



iBoogert et alj (|2002|) illustrates how complex the environment ridge of material, but in front of that ridge, two foreground lay- 
can become. Elias 29, an embedded YSO, is located in a dense ers were identified using the emission from a dozen different 
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Fig. 13. The spectra of C'^O 3-2 and HCO+ 4-3 of the LFAM 26 field. The intensity scale is -1 to 6 K in Tmb for C'^O and -1 to 
3 K for HCO^. OphE-MM3 is located at ARA = +5" and ADec = -40" in this field and has no associated HCO^ emission, except 
from the extended cloud. 
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Fig. 15. The spatial distribution of the integrated gaussians fitted to the C'^O 3-2 spectra for each of the four layers seen in the field 
of LFAM 26. 

Table 5. The number and properties of the (foreground) layers for the fields of LFAM 26, IRS 44 and WL 19 in the Ophiuchus ridge 
and IRS 37 and IRS 4, located near Oph-B cores, see Fig. |7] 



Field 


Position 


Nr. comp 


Vlsr (kms-') 


C'^'O j'TMBdV(Kkms-^ 


Ay(kms"') 


LFAM 26 


0,0 


3 


3.5,4.4,5.0 


1.8, 0.9, 1.1 


0.8,0.7,0.8 




-20,+35 


4 


2.7,3.5,4.4,5.0 


0.6, 1.5,0.8, 1.1 


0.9, 0.7, 0.7, 0.7 




+35, -22 


3 


3.5,4.4,5.0 


1.3,0.4, 0.5, 1.2, 0.7,0.5 




IRS 44 


0,0" 


3 


2.6, 3.5, 4.5 


0.5,2.3,2.3 


1.9,1.1,0.8 




+20,+20 (IRS 46) 


2 


3.5,4.3 


1.9, 1.3 


0.6,0.5 




-50,20 


2 


3.3, 4.0 


1.1,2.5 


0.4, 1.1 


IRS 48 


0,0 


3 


2.8,3.5,4.1 


0.6, 1.7,2.6 


0.6, 0.4, 0.5 




-36,+36 


3 


2.9,3.5,4.1 


0.8,2.1,4.3 


0.4, 0.4, 0.5 




+38, 6 


3 


2.8,3.5,4.1 


0.4, 1.1,3.5 


0.3, 0.5, 0.5 


WL 19 


0,0 


2 


3.3,4.6 


2.7,4.7 


0.8,0.8 




-6, 36 


2 


3.4, 4.6 


2.8,5.6, 1.1,0.8 






-20,-38 


2 


3.2,4.2 


1.9,3.0 


0.7, 1.0 




-50, 50 


4 


2.6,3.5,4.3,5.1 


1.4,3.3,3.5,3.2 


0.5,0.6, 0.7,0.7 


IRS 37 


0,0 


2 


3.0, 4.4 


2.3,5.6 


0.7, 1.1 




+23, +7 (WL3) 


2 


3.2,4.4 


2.7,5.4 


0.8, 1.0 




-36, 21 


3 


3.0,3.6,4.2 


2.4, 0.8, 2.8 


0.9, 0.2,0.8 



"An outflow is detected at 0,0 in HCO^ 4—3 which is probably the component seen at 2.6 km s 



molecules. Continuum emission and ice absorption originate in 4-3 profiles are mostly single peaked. However, C"*0 3-2 often 
all these layers. For most of our sources, however, the situation shows a more complex profile with multiple peaks, 
appears to be less complex. Fig. |6] and |7] show that the HCO^ 
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Fig. 12. The distribution of concentration of SCUBA (gray) and 
HCO^ 4-3 (striped). Only sources with spatially resolved pro- 
files have been taken into account. Note that the total number of 
sources is greater for SCUBA than for HCO^. 
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Fig. 14. Individual spectra at offset positions from LFAM 26 
for HCO+ (0",0") and C'^O (0,0), (-50,50), (10,50) and (25,- 
40). The three main components at 3.5, 4.4 and 5.0 km s"' are 
shown with solid lines. The component at 2.7 km s ' identified 
by Boog ert et al.i ( f2002l) is shown with a dashed line and is gener- 
ally weaker or absent. The main component of the HCO^ emis- 
sion is located between the components at 4.2 and 4.9 km s"^ 
and peaks at 4.6 km s ' . 



LFAM 26, located at the top of the Ophiuchus ridge only a 
few arcmin north-east of Elias 29, is an especially complex case. 
Fig.[T3]shows the spectral map in 17" spatial bins for HCO^ 4-3 
and C'**0 3-2 maps. It is clearly seen that besides the integrated 
emission (Fig. [3]), the line profiles of HCO^ and C''^0 both vary 
over the entire map. The HCO^ spectra show mostly single gaus- 
sians. At the position of LFAM 26, the HCO^ profile has a small 
outflow signature on the blue side of the spectrum. The C'^O 
emission requires at least four distinct components (Fig. [T4l i at 



2.7, 3.5, 4.4 and 5.0 km s The three c omponents at 2.7, 3.5 
and 5.0 km s"' are the same as found in 'Boogert et al." (20025 
at the position of Elias 29. Most lines observed in Boogert et al] 
(|2002|) peak short-wards of 5.0 km s ' component, which could 
be indicative of the 4.4 km s ' component . However, the differ- 
ence between these two components is very pronounced in the 
spectra in Fig. [13] and [14] at the position of LFAM 26 so that 
we treat them as separate. HCO^ is aligned with the 4.4 km s"' 
component. 

Using the molecu lar excitation program RADEX 
(Ivan der Tak et all l2007l) . one can calculate the expected 
average contribution from cloud material to the HCO^ 4-3 
emission. Using a HCO^ abundance of 10 with respect to H2, 
a density of lO'* cm~^, a cloud temperature of 15 K, a width 
of 1 km s ' and a H2 column density of 10^^ cm"^ (Ay ~ 
10), a contribution of ~ 0.1 K km s ' is found. Only at higher 
column (~10^-' cm"^) or volume densities (~lxlO'' cm"-') are 
contributions found from the cloud on the order of 1 K km s 
comparable to the observed 4.4 km s"' emission at off-positions 
such as (-50, 50). Th i s is in agreement with the numbers given 
in Table 1 of lEvansI (1 19991) . Thus, LFAM 26 is embedded in a 
high-density ridge. 

Table|5]shows the results of gaussian fitting at other positions 
where multiple peaks in the C'^O 3-2 profiles were found. In 
the WL 19 field, south of the Elias 29 position, the same four 
components are again detected, although the 5.0 km s ' layer is 
only seen north of WL 19. In the IRS 44 field along the same 
ridge to the south-west, the three layers at 2.6, 3.5 and 4.5 km 

are found. 

Three layers are found in the Oph-B2 core, which contains 
the IRS 48 and IRS 37 sources. The first layer is at 2.8 km s ', 
the second at 3.5 km s ' and the third at 4.2 km s '. Outside 
these two regions, spectra of C"^0 can be fitted with a single 
gaussian, almost always at ~ 3.6 km s"'. The spatial variation 
in this layer seems to be minor compared with the other compo- 
nents and it is thus the most likely candidate for the foreground 
component with an extinction of Ay ~5.5 mag, as suggested by 
iDickman & HerbstI ([1990). IRS 63, located far from the L 1688 
core is at 2.8 km s"\ the same velocity as one of the layers in 
both Oph-B2 and Ophiuchus ridge region. 

5.3. Gas column density 

The total column density, A^h,. toward each source was calcu- 
lated from the integrated intensity of the C"*0 3-2 line,JrMB 
dV, using the following formula : 

[TMBdV 

NiH2) = X X 0.75 X 10'^ i_,r-33.2/r) (3) 

with X, the abundance of H2 with respect to C'^O, assumed to 
be 5.6x10'' ( Wilson & Rood 1994) and T, the temperature of the 
cloud material, assumed to be 30 K. The temperature is chosen 
to be higher than that of the surrounding cloud material, ~15 K, 
and reflects the warmer gas associated with the YSO where CO 
is not frozen out. This formula assumes the molecular excitation 
to be in LTE and isothermal. The C'^O 3-2 line is assumed to 
be completely optically thin. The resulting column densities can 
be found in Tabled] column 3. They range from 3.7x10^^ cm"^ 
for IRS 63 and the Oph F core to 3 1 .6 x 10^^ cm'^ in the Oph A 
core. These correspond to visual extinctions Ay of 40-400 mag_^, 
assuming Na/Ay =1.8x10^' cm"^ mag"' (Rachford et al. 20()2) 
and A^H ~ 2A^(H2) (all hydrogen in molecula r form). This range 
of column densities is similar to that found bv lMotte et an(ll998l) 
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and IStanke et alj ( l2006l) using 1.3 mm continuum imaging, as 
well as the column densities found from the 850 yum dust, see 
§5.4. 

5.4. Dust 

Table|4]includes the properties that can be derived from the con- 
tinuum emission. The dust column density was calculated with 
the following equation. 



A^(H2) 



> 850 



n ju mn KS5Q Bs5oiTd) 



(4) 



where S sso is the flux in the central beam, Q is the main beam 
solid angle of 15", fi = 2.33 the mean molecular weight, niu the 
mass of atomic hydrogen, /cgso the dust opacity per gram of gas 
and dust at 850 fim and Bg5o(r) the Planck function at 850 fim 
for a temperatu re T, assumed to be 30 K. A Kp p - 0.01 cm^ gr"' 
is adopted from lOssenkopf & Hennin j (1 19941) for dust with thin 
ice mantles (Table 1, column 4). 

Two assumptions affect the comparison between the column 
densities derived from C"^0 and those from dust continuum. 
C'^O can be frozen out onto the dust grains below ~30 K 
at sufficiently high densities, so that C"^0 column densities 
are expected to be lower than those obtained from dust data. 
Second, all dust emission is assumed to originate in envelope or 
cloud material, whereas for several sources the disk contributes 
significantly to the continuum. These combined effects likely 
eplain the lack of correlation seen in Fig{T6](/e/f) 

If one assumes that all emission detected at 850 /^m is 
associated with an isothermal envelo pe, a source mass ca n be 
determined following Equation 4 from lShirlev et alj (l2000l) : 



S 850 D 

By(Td) Ky 



= 3.69 X 10" 



x585o(Jy)D2(pc) X (e 



16.9/Td 



Me 

- 1) 



(5) 



Since the 3cr sensitivity of the COMPLETE map is 80 mJy 
within a 15" beam, the lower limit to the detectable mass is 0.04 
Mq. Such a low mass implies either a source that has nearly shed 
its envelope and is close to the pre-mai n sequence phase or an 
intrinsically very low luminosity object jPunham et al.ll2008h . 

Fig-dHshows the relation between the different column den- 
sities, envelope masses and the presence of dense gas as traced 
by HCO^. These figures clearly show that at higher envelope 
masses, the amount of dense gas and the total envelope mass cor- 
relate. The third figure illustrates that nearly all sources with cen- 
trally concentrated HCO^ emission have larger envelope masses. 
Above an envelope mass of 0.1 Mq, only a single source that is 
classified as confused in § 6 (indicated with diamonds) is seen 
with HCO+ > 2 K km/s. This is CRBR 2422.8-3423, which is 
located in the line of sight with the envelope of IRS 43. The high 
mass found for CRBR 2422.8-3423 can be fully attributed to this 
envelope. 

5.5. Effect of reddening on aiR 

Foreground material reddens the YSO and steepens the spec- 
tral slope Q'2-24/jin- An estimate can be made of the foreground 
column density based on our C'**© maps, by assuming that the 
average column or Ay at 45" offset from the source position in 
all directions originates entirely from this cloud. Although it is 



not known a priori whether this material is foreground, back- 
ground, or a different YSO, we assume here that it is all fore- 
ground (see also §5.2) which provides an upper limit to the red- 
dening. Note that the column densities at these positions were 
calculated with a temperature of 15 K, a temperature associated 
with the surrounding cloud, instead of 30 K, a temperature as- 
sociated with protostellar envelopes. For typical du st opacities 
and the Ay-5.5 mag layer as suggested by Dick man & Herbsl 
(1990), a2-24/jm is only increased by 0.18 (see also lCrapsi et al. 
(i2008i) ). However, for higher Ay of ~ 40 mag, as commonly 
found in the Ophiuchus ridge, this correction can be as high as 
1.1, enough to change the classification of a Class II T Tauri 
star with disk to an embedded Class I source. Table |4] includes 
the (maximum) change in a2-24;im due t o foreground reddening 
for each source derived by this method. ^ Chapman et al.l (l2008l) 
probe the Ak in the northern part of Ophiuchus, including the 
northern part of the Oph-A area, where Ay values of 40 mag. 
and higher are found. In the Oph-A, Oph-B2 and Ridge regions. 
Ay of this magnitude were found from C'^O column densities, 
reinforcing the conclusion that a strong extinction affects these 
regions on scales of only a few arcminutes, a resolution often not 
reached by extinction studies dCambresvll 1 9991) . Although these 
numbers have large uncertainties, they do illustrate the point that 
this potential extreme foreground reddening can lead to incorrect 
classifications for a significant number of sources in Ophiuchus. 

An alternative is to use the velocity resolved C'^'O spectrum 
at the source position but subtract the layer that includes the 
source. For example, LFAM 26 (Aq'ir=4.4), resides in a dense 
ridge (see § 5.2). Subtracting the gaussian profile of this layer 
at ~4.4 km s"' gives an integrated intensity of 2.9 K km s"', 
which corresponds to a column density of 10.7x10^^ cm"^, an 
Ay of ~60 mag and Aq-ir = 3.2. If one assumes that the layers 
are arranged in distance in order of increasing Vlsr, only the 
front layer at 3.5 km s ' will redden and a column density of 
6.6x10^^ cm"^, an Ay of ~40 and Aci'ir=1.9 are found. Thus, 
AaiR is likely to be in between the maximum value listed in 
Table 4 and Aq'ir/2, although exceptional cases, such as OphE- 
MM3, an edge-on disk in front of the cloud, do exist. 



5.6. SEDs;Lboi and Tboi 

The IR and submillimeter fluxes have also been used to calculate 
the bolometric luminosity, Lboi, and temperature, Tboi, of each 
source. These numbers were calculated using the prismodial or 
the midpoint methods, with the inclusion of the 2MASS, IRAC, 
MIPS-24, MIPS-70, SHARC-II, SCUBA and 1.3 mm fluxes. 
Only fluxes with aS /N > 5 were included. For a more thorough 
discussion on the useage of these two methods, see Dunham et 
al. (in prep.) and .Enoch et al., (2008). 

Table |6] compares the values found for Tboi and Lboi in 
our work with those of lEvans et al.l (|2009|) for several sou rces. 
The differences stem from the fact that E vans et al] (|2009|) use 
CSO-Bolocam 1 . 1 millimeter data in a 30" beam, together with 
2MASS, IRAC, and MIPS (24, 70 and 160 /im) fluxes. No 
SHARC-II or SCUBA fluxes were included, nor were sources 
with no detection at 1.1 mm. As can be seen from the results 
for Elias 29 and GSS 30 IRSl, the inclusion of Bolocam and 
MIPS-160 produces a higher Lboi of up to a factor of 6. In ad- 
dition, Jboi is consistently lower for the brighter sources. This 
can be largely attributed to the contributions from the surround- 
ing cloud, which add significantly to the MIPS-160 and Bolocam 
fluxes. Although cloud emission is present in the IRAM-30m 1.3 
mm and SCUBA-850 fim observations used in our studies, the 
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Fig. 16. Le/f : the gas column density (in 10^^ cm"^) derived from C'^O 3-2 vs. the envelope mass (in 10"^ Mq), derived from the 
SCUBA maps; middle : the column density vs. the integrated intensity of HCO^ 4-3 in K km s ' and right : the integrated intensity 
of HCO^ 4-3 vs. the envelope mass . In all figures, sources with a centrally condensed HCO^ emission are shown with stars, sources 
with an upper limit on HCO^ 4-3 with triangles, and sources with HCO^ emission not peaking on the source with diamonds. 



higher resolution of 15" limits such contributions. The exclusion 
of MIPS- 160 in our work causes the Lboi to be underestimated 
and the Tboi to be overestimated, due to the lack of a point at far- 
IR wavelengths near the peak of the SED. Isolated, unresolved 
sources such as IRS 42 show little to no differences between the 
two studies, confirming the influence of the environment. Future 
high resolution (<10") observations with the PACS instrument 
on Herschel covering the 60 to 600 yum window will be able to 
fully constrain the far-IR emission and SED. 

Table 6. C omparison between Lboi and Tboi as calculated in 
lEvans et aTl (12009.) (indicated with Evans) and this work (indi- 
cated with van K.). 



Source 








Tboi 




Evans 


van K. 


Evans 


van K. 








(K) 


(K) 


GSS 30 


8.7 


3.3 


150 


123 


CRBR 2339.1 


<0.065 


0.055 


<370 


426 


WL2 


0.12 


0.09 


428 


573 


LFAM 26 


0.15 


0.044 


110 


238 


WL 17 


0.6 


0.67 


310 


323 


Elias 29 


17.9 


2.5 


257 


424 


IRS 42 


1.2 


1.4 


600 


540 


CRBR 2422.8 


0.19 


0.44 


<300 


157 


IRS 43 


3.8 


1.0 


160 


134 


IRS 44 


15 


1.1 


110 


140 



5.7. Ice and silicate absorption 

In the online appendix, the IRS spectra for 28 of the sources in 
our sample are plotted. 17 sources show ice absorption bands 
at 15.2 fim due to CO2 ice (see Table |7]), and 11 do not (WL 
3, LFAM 26, VSSG 1, CRBR 2317.3, CRBR 2315.8, IRS 48, 
CRBR 2339.1, SR 21, Hai-o 1-4, IRS 46 and SR 9). A few of the 
known edge-on disks, such as CRBR 2422.8-3423.2, show ice 
absorptions. 

ICrapsi et alJ ([2008') modelled the spectral features at 3 fim 
(H2O ice), 10 fim (silicate) and 15.2 //m (CO2 ice) for a large 
grid of models, including embedded sources and T Tauri disks, 
both seen at a wide variety of angles. The ice absorptions at 3 and 
15.2 fj.m were found to be more prevalent for embedded sources. 
However, edge-on disks can also show such absorptions in al- 
most equal strength. The relation of the 3 //m ice absorptions 



with the envelop e mass of bo t h emb edded and disk sources is 
given in Fig. 4 o f ICrapsi et alJ (l2008h and a similar relation was 
found for the 15.2 fj.m band. Cold foreground clouds can also 
contribute significantly. 

Using the determination of the column densities from the 
C'^O and SCUBA maps in § 5.3 and 5.4, Eq. 3 and 4, the ice 
absorption can be compared to the amount of reddening as in- 
dicated by AffiR. Although ice absorptions are more commonly 
found in sources with a higher reddening, it is concluded that the 
YSOs must be characterized on a source by source basis to locate 
the origin of the ices. The disk source C RBR 2422.8 has deep ice 
absorptions dPontoppidan et al] l2005h . most hkely originating 
the large column of foreground material (Aa = 0.9). However, 
LFAM 26, which has an even higher reddening (Aa - 1.7), 
shows no sign of ice absorption. Of course, it is possible LFAM 
26 is located in front of the cloud and as such the presence or ab- 
sence of ice absorption could serve as a useful diagnostic of the 
geometry of the YSO-cloud system. The silicate feature around 
10 jum is detected for fewer sources than the 15.2 CO2 ice ab- 
sorption. Although LFAM 26, WL17, Elias 29, IRS 42, WL 3, 
IRS 54, CRBR 2315.8-1700, CRBR 2339.1-2032 and CRBR 
2422.8-3423 all have silicate absorption, many other sources 
were not observed at the wavelengths of the silicate feature. 

6. Classification 

6.1. Physical classification 

A new classificat ion based on physical p a rameters has gradually 
been introduced (IWhitnev et al.ll2003by. [Robitaille et alj|2006l) 
but is not yet as commonly in use as the traditional Class system. 
This new classification identifies the evolutionary stage using the 
physical characteristics as opposed to the observed characteris- 
tics. The different evolutionary stages are determined from the 
ratios between Mdisk, A^enveiope and Mstar. The total circumstellar 
mass, Mcircum, is defined as Mdisk+A^enveiope- The evolutionary 
stages of this classification are: 

- Stage 0, deeply embedded sources with Mdisk/Mgnveiope << 1 

and Mcircura/Mstai- ~ 1 

- Stage 1, embedded sources with 0.1 < Mdisk/A^enveiope < 2 

and Mcircura < Mstai- 

- Stage 2, classical T Tauri stars with gas-rich disks (Menveiope 
= and Mdisk/Mstar « 1 ) 

- Stage 3, Pre-main sequence stars with tenuous disks. 

Stage sources are equivalent to Class sources and can be 
identified by their submm characteristics as originally put for- 
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Table 7. The classification of the source sample with the recipe 
outlined in the text. 



Source 


Code" 


Stage 


Embedded sources 


GSS 30 


Shcs E/f8 PhsI 


1 


LFAM 1 


S„8 Nc E„8 P„ Og 


1 


VLA 1623 


ShC8 E/yg P;/8 




WL 12 


Shc8 E/ycg PhcsI 


1 


LFAM 26 


Shcs E/^cs Phs Oc 


1 


WL 17 


Shs Nc P/fsI 


1(T) 


Elias 29 


S/yCS E/y8 P/igl 




IRS 37 


Shc Wg Emcs PhcsI 




WL3 


Shc Wg Ehcs Ph 




WL6 


Shc Ngl 


1(T) 


IRS 43 


Shcs Ehs PhsI 




IRS 44 


Shcs Ehcs Phcs 




Elias 32 


Shcs Ehs PhsI 




Elias 33 


Shcs Ehs PsI 




IRS 54 


Shc Wg Efjcg PhcsI 


1(T) 


IRAS 16285-2355 


Shs Wc Emcs Phcs 




IRS 63 


Shs Wc Emcs PhcsI 


1(T) 


Confused Sources 


CRBR 2324.1-1619 


Shcs Ehcs Ohcs 


C 


GY51 


Ss Nhc Eg Og 


C 


WL19 


Scs N// Efjcs Ohcs I 


c 


IRS 42 


Shcs Ehs OhsI 


c 


GY 256 


Shc Ng 


c 


Disks 


SSTc2d J162527.6-243648 


Nhcs 


2 


GSS 26 


Ss Nhc Ps 


2 


CRBR 2315.8-1700 


Shs Nc E„ Pg Oh 


2 


CRBR 2317.3-1925 


Sh Wg Nc 


2 


VSSG 1 


Wg N„c Ps 


2 


CRBR 2339.1-2032 


Sh Wg Nc Pg 


2 


WL2 


Whs Nc Ehs Ohs 


2 


GY 224 


Wg Nhc PsI 


2 


WL 20S 


Scs Nh PsI 


2 


IRS 48 


Wcs N« Ec Pg 


2 


GY312 


Nhcs 


2 


IRS 51 


Shs Wc E^cs Pg OhcI 


2 


SSTc2d J162741.6-244645 


WgN«c 


2 


SSTc2d J162748.2-244225 


Nhcs 


2 


SSTc2d J162857.9-244055 


Ss N„c 


2 


Haro 1-4 


Nhcs 


2 


DoAR 25 


Sg N//C Pg 


2 


OphE MM3 


Shc Wg Ehcs Ohs 


2 


SR21 


Wg N„c Ps 


2 


CRBR 2422.8-3423 


Shs Nc Eg Ogl 


2 


IRS 46 


Ss Wc Nh Ehcs Ohcs 


2 


SR9 


Nhcs 


2 


2Mass 16282 


Wg Nhc 


2 



"The coding in column 2 is as follows: 

S, W and N determine if a line or continuum is Strong, Weak or Not 
detected/observed. 

E, P and O determine if a source is Extended, Peaking or Offset peaking. 
The subscripts H, C and 8 refer to the HCO+, C'**0 and 850 ijm. 
I is detected CO2 ice absorption. Column 3 lists the classification. 1 is a 
Stage 1 embedded YSO, 2 is a disk and C is Confused. 



ward by lAndre et al.l (Il993h . Stage 2 and Stage 3 sources, corre- 
sponding to the classical T Tauri stars with disks and pre-main 
sequence stars with (gas-poor) tenuous disks, are identified using 
the IR excess combined with optical/near-lR spectroscopic prop- 
erties. The amount of IR excess and the wavelength where this 
excess starts are often used to distinguish Stage 2 and 3 sources. 



The Stage 1 sources are most problematic to uniquely iden- 
tify based on observational characteristics. Both IWhitnev et all 
(2003a) and Crapsi et al. (2008) show that the traditional iden- 
tification methods of using the IR spectral slope ffiR and Tboi 
are insufficient to distinguish the edge-on disks from embedded 
sources, as well as missing evolved, face-on embedded YSOs 
that are (mis)classified as Class II. One of the main parameters 
to unambiguously identify Stage 1 sources is the Mdisk/A^enveiope 
ratio as determined from milli meter continuum o bservations in a 
large and small (~ 2") beam dCrapsi et al.ll2008 ). However, the 
interferometric observations needed to unambiguously constrain 
the disk masses are still time-consuming for large samples of 
sources. 

6.2. Identifying embedded stage 1 sources w/fh moiecuiar 
emission 

To uniquely identify embedded sources and distinguish them 
from edge-on disks, prestellar cores and possible background 
AGB stars and galaxies without having to resort to interfer- 
ometers, we propose to use the single-dish molecular emission. 
Embedded sources will be bright and centrally peaked in HCO^ 
4-3 and C'**0 3-2. The HCO+ Hne maps will detect all but the 
most tenuous envelopes. The concentration parameter shows that 
HCO^ 4-3 is much more accurate in identifying protostellar en- 
velopes than continuum dust emission. The HCO^ maps also 
spatially resolve larger envelopes (e.g., the VLA 1623 and GSS 
30 envelopes). The nearby environment, including foreground 
layers, is in turn characterized by the C"*0 3-2 (see § 4 and 5). 

To distinguish embedded sources from isolated edge-on 
dis ks, a good limit on the HCO^ 4-3 emission is Tmb > 0.4 
K. iThi et alj ( |2004|) observed four well-studied gas-rich disk 
sources in a large range of molecular transitions, including 
HCO^ 4-3, and found Tmb ~ 0.1 K, four times below our 
adopted limit. Although disks around embedded sources could 
be somewhat larger and more massive, it is assumed that the 
contribution from even the largest embedded disks does not ex- 
ceed our adopted limit. This value is naturally small because of 
the severe dilution of emission originating in a disk within the 
15 " JCMT beam. For example, a 20 K optically thick line from 
a 2" disk would be diluted to 0.35 K. Indeed, for two embedded 
sources f or which interferom etric HCO^ data exist (IRS 63 and 
Elias 29, iLommen etakl (I2OO8)), the HCO+ contribution of the 
disk does not exceed this limit (see §6.3). 

The following definitions were therefore used to classify the 
sample: 

- Stage 1 

- Extended HCO^ emission, peaking on source, with 
Chco* > 0.7 and J TuBdV > 0.4 K km s^' and/or 

- an HCO+ profile that is not extended, but J TMBdV > 0.4 
K km s with SCUBA both extended and Cggo > 0.7. 

- Ice absorptions are usually prominent. 

- Stage 2 

- Absence of HCO+, SCUBA and C"^0 down to the rms 
limits at the IR position or 

- no variation of HCO+, SCUBA and C'**0 on scales of 
30" (C < 0.4) or 

- no C'**0 3-2 or SCUBA is seen at scales larger than the 
central beam of 15" with HCO^ 4-3 not extended and < 
0.4 K km s"' or 

- aiR < -0.5. 

- Confused 
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- HCO+, SCUBA or C'**0 peaking more than 20" away 
from the IR position 

In the rare case that a n embedded source is viewed face-on 
(into the outflow cone, see lWhitnev et alj (l2003bh ). the above re- 
strictions will identify such a source as stage 1, having am < 
-0.3, but with extended HCO^ emission, peaking on source 
(/ rMBdV>0.4Kkms-i). 

Table |2] shows the results for this classification. In Appendix 
C the sources and their classifications are discussed on a source 
by source basis. Sixteen Class I sources were identified as em- 
bedded Stage 1 sources with envelopes varying in mass and size. 
The sources GSS 30 IRS 1 and LFAM 1 appear to be embed- 
ded within an approximately spherical envelope encompassing 
both sources. The sources IRS 37 and WL 3 are embedded in a 
highly non-spherical envelope covering both sources. In addition 
to these 16 sources, VLA 1623 is classified as embedded, but 
considered a Stage source due to its high sub- mm flux relative 
to the total luminosity. From comparison with |j0rgensen et al.l 
(I2008h (see § 6.5), we can conclude that all embedded sources 
were found down to a mass limit (envelope+disk) of 0.04 Mq. 
Deep searches for even lower luminosity sourc es have found 
only t wo likely candidate sources within our field jPunham et al.l 
|2008|) . Deep (interferometric) sub-mm observations are needed 
to identify if such sources are truly embedded. 

Five sources were found to be confused (C in Table |7]i. 
The exact nature of these sources (obscured, disk or other) 
cannot be determined without interferometric observations. 
However, molecular emission rules out a protostellar envelope. 
Observations of C'^O show that cloud material is located in the 
line of sight towards these sources. They will be further dis- 
cussed in §6.4. 

Twenty-three sources were found to be Stage 2 sources (2 in 
in Table |2ll. Of this sample, sources with cir > and with no 
correction from extinction (see Table 4) are likely to be edge- 
on or close to edge-on disks. Sources that do have a potential 
correction to ajR, such as CRBR 2422.8 or IRS 48, do not have 
to be edge-on. All sources that were initially known to be disk 
sources are identified as Stage 2 sources using this method. If 
one considers that disk sources with aju >-0.3, when corrected 
for the maximum Aq-ir), are edge-on, only 7 out of 23 Stage 
2 sources of our sample are expected to be edge-on. The total 
amount of Class II sources found in Ophiuchus by c2d is 176, 
and although this includes L 1689, the bulk of these disks exis t 
close to or are part of the L 1688 region (lEvans et al.ll2009l) . 
Thus, our number of edge-on disks appears to be only a small 
fraction of the total disk sample. 

No face-on embedded sources were found that have q-ir < 
-0.3, but do show a HCO^ spectrum that is both peaking on 
source and stronger than 0.4 K km s This very low fraction is 
not surprising, since such sources have strong restrictions on the 
line of sight. 

Using the result s from Spitzer jEvans et aL I (I2007h . 
|j0rgensen etaTI (l2006h . lEvans et al.1 (l2009h . Allen et al. in prep) 
a comparison can be made between the classification using q-ir 
and the new method. It is found that ~50% (11 out of 22) of the 
sources classified as Class I are Stage 2 disks and not embedded. 
Note that this includes the known disks OphE MM3 and CRBR 
2422.8. However, of the Flat spectrum sources, ~5Q% (6 out of 
13) appear to be embedded, while the other half is either con- 
fused or classified as a Stage 2 disk. In the end, the sample of 
embedded sources found by c2d changes from 22 out of 45 to 
17 out of 45. 



Table 8. Comparison bet ween the classification of YSOs using 
aiR bv lEvans et al.l (l2009l) and our method using molecular emis- 
sion. 





Class I 


Flat Spectrum 


Class II 


Total 


Stage 1 


11 


6 





17 


Confused" 


1 


3 


1 


5 


Stage 2 


10 


4 


9 


23 


Total 


22 


13 


10 


45 



"Note that although it is possible for Flat Spectrum sources to be em- 
bedded, confused sources are ruled out to be embedded. 



6.3. Late stage 1 sources 

Of the sixteen embedded Stage 1 sources, four sources (WL 6, 
WL 17, IRS 54 and IRS 63) were found to have only marginal 
envelopes left, with C'**0 and 850 fim emission suggesting that 
at mos t a few xO.OlM© remains in the envelope. iLommen et al.l 
(12008 ) observed IRS 63 using the Sub-Millimeter Array (SMA) 
and found that 50 % of the continuum emission originates within 
a disk no bigger than 200 AU in size with a mass of 0.055 Mq 
and the remainder in a remnant envelope with a mass of 0.058 
Mq. HCO^ 3-2 was also detected with the SMA, originating 
mostly in the disk. When convolved with the JCMT beam, the 
observed brightness is 0.28 K km s The HCO^ 4-3 single- 
dish intensity as observed with HARP-B is 0.75 K km s 3 
times as bright. For typical disk excitation conditions, HCO^ 4- 
3 is expected to be equal or weaker compared to the 3-2. It can 
thus be concluded from the much higher observed intensity in 
the 4-3 line, that at least half of the 4-3 emission within the 
JCMT beam likely originates in the protostellar envelope and 
would not be seen by the SMA. 

In all four sources, little to no extended C'**0 is seen. The 
distributions and strengths in the HCO^ and C'^O lines of the 
other three sources are similar to that of IRS 63, but the SCUBA 
image of IRS 63 is a factor of 4 brighter For IRS 54, the enve- 
lope is spatially resolved in HCO+, C"^0 and SCUBA. 

Although rare, it is concluded that these sources are Stage 
1 sources in transition to a Stage 2 source (marked with 1(T) 
in Table |7]i. They have accreted or dispersed almost their entire 
envelope. Continuum emission in the sub-mm is likely to contain 
a large contribution from the protostellar disk. 

6.4. Confused sources 

For the five confused sources, CRBR 2324.1-1619, GY 51, WL 
19, IRS 42 and GY 256, the exact nature cannot be determined 
using HCO^, except that such sources are ruled out as embed- 
ded YSOs. None of them have the characteristic peak of HCO^ 
that most Stage 1 embedded YSOs have. However, at least two 
out of the three tracers, HCO+, C'**0 and SCUBA, are detected 
so a few possibilities remain. The first is an (edge-on) disk, in 
front of cloud material. This would be similar to the case of 
OphE-MM3, which was confirmed to be an edge-on disk in 
front of a dense core (' Brandner et al.ll2000 ^. Without the near- 
IR imaging, OphE-MM3 would have been classified as con- 
fused. A second possibility is that they are background sources. 
The most likely options are then a T Tauri star with disk be- 
hind the cloud (possibly edge-on), or an AGB star Background 
main sequence stars are identified within the c2d delivery doc- 
ument (Evans et al. 2007) and are thus highly unlikely to be 
included in our sample. Chances of background AGB stars or 
galaxies being aligned with the cloud are small, but according to 
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lJ0rgensen et alj (l2008l) . a single background source can be ex- 
pected. A final possibility is that the sources are very late Stage 
1 embedded YSOs, very close to the Stage 2 phase. An upper 
limit of only 0.02 is found for their envelope mass. 

6.5. Comparison to other methods 

Fig. [17] shows the effectiveness of the classification of sources 
using the method above, as compared to methods using q-ir and 
Tboi. A key parameter is the strength of the HCO^ integrated 
intensity. The following limits were adopted between embedded 
and disk sources: ckir > 0, Tboi < 650 K and HCO+ > 0.4 K km 
s-i. 

The advantages of the method using molecular emission are 
immediately apparent. Although the classical methods of using 
Tboi and am are able to identify embedded sources, both meth- 
ods also incorrectly identify a number of Stage 2 disks as embed- 
ded, especially for Lboi < 1 Lq. These are the edge-on disks. If 
one uses the limits of the integrated intensity for HCO^, isolated 
disk sources are easily identified. With the additional restriction 
of having a peak within the HCO^ map, confused sources or 
sources in front of the cloud can be readily identified. 

The use of ice absorption, as studied bv lCrapsi et al.l (l2008h . 
is limited by similar constraints. Most embedded sources, with 
the exception of LFAM 26 and WL 3, show CO2 ice absorption 
at 15.2 fj.m. However, 4 disk sources IRS 51, WL 19, WL 20S 
and CRBR 2422.8-3423 also show strong ice absorptions. These 
ice absorptions are caused both by the foreground layers as dis- 
cussed in §5.2 as well as the material in the disk itself, if viewed 
edge-on (as is the case for CRBR 2422.8). Thus, the presence of 
ice absorption cannot be used to unambiguously identify embed- 
ded sources from other sources. Even if the origin of the ice ab- 
sorption can be attributed to a protostellar envelope, its strength 
does not seem to be an indication of evolution, e.g. IRS 63, a 
small envelope, has a deeper ice absorption than the large enve- 
lop eofGSS30. 

|j0rpensen et al.l (l2008l) published a list of 'candidate' embed- 
ded YSOs based on two criteria. First the colors of sources using 
the IRAC and MIPS results ([3.6] - [4.5] and [8] - [24]). Second, 
the proximity of MIPS sources to SCUBA cores. For these cores, 
the concentration is an important parameter. Although most can- 
didate embedded objects are identified by both criteria, several 
sources are included in their list based on only a single one. 
Comparison between their list (Table 1 in J0rgensenet al. (2008|) 
limited to L 1688, further referenced as JJI) and the list in this 
paper yields the following results: 

- Four embedded YSOs with SCUBA fluxes below the cut-off 
adopted by JJI (< 0.15 Jy beam the very low envelope 
masses) are absent from the JJI list. Two of these sources 
(WL6 and IRS 54) are classified by us as late Stage 1 sources 
with little to no envelope left. The other two embedded YSOs 
not included in JJI are IRS 37 and WL 3. 

- Four sources are included in the JJI list that have been clas- 
sified as Stage 2 disks by our method. These are GSS 26, 
CRBR 2315.8, CRBR 2339.1-2032 and IRS 51. AH four 
sources have associated SCUBA cores and MIPS detections. 
See the appendix for the classification reasoning of each of 
these sources. 

7. Conclusions 

A sample of young stellar objects in L 1688 was analyzed using 
gas mapping obtained with the new HARP-B heterodyne array 



receiver in the HCO^ 4-3 and C'^O 3-2 lines. Complementary 
dust maps were obtained from the COMPLETE project, as ob- 
served by JCMT-SCUBA, and with SHARC-II on the CSO. The 
original sample consisted of 45 sources, mostly classified as em- 
bedded YSOs or flat-spectrum sources using their spectral slope 
aiR in previous work. Of this sample a few sources were recently 
discovered to be edge-on disks. As a control sample, 4 known 
disk sources in L 1688 were included. The observations were 
supplemented by single-pixel observations from APEX, Spitzer- 
IRS spectroscopy and continuum photometry ranging from 1 jum 
to 1.3 mm, using a variety of space-based and ground-based ob- 
servatories. 

The main conclusions are: 

- The concentration of the dense gas, as traced by the HCO^ 
4-3 line mapping, provides an excellent tool to characterize 
the dense gas in the inner regions of protostellar envelopes. 
Material in the cold outer envelopes, (edge-on) disks, prestel- 
lar cores or cloud material does not emit strongly in HCO^ 
4-3. 

- Most envelopes in L 1688 have low masses, ranging from 
0.05 to 0.5 Mq. The main accretion phase onto the star has 
already taken place. The only exception is the Stage source 
VLA 1623 which contains nearly 1 Mq. 

- A new classification tool based on molecular emission to 
uniquely identify embedded sources is proposed. Of the total 
sample, 17 sources were found to be embedded, 5 sources 
were confused and 23 sources were identified as (edge-on) 
disks. Due to the foreground layers found through the C'^O, 
it is concluded that a significant fraction of the sources pre- 
viously classified as Class I or Flat Spectrum are disks and 
have been strongly reddened. Combined with the analysis in 
|j0rgensen et al. (200§) all embedded sources in L 1688 were 
included and characterized down to a mass limit of 0.04 M© 
(envelope-i-disk). 

- Four embedded YSOs are found that are in transition from 
the Stage 1 embedded phase to the Stage 2 T Tauri phase. 
These sources, IRS 63, IRS 54, WL 6 and WL 17, have litde 
envelope left (Menv ~ O.O5M0). 

- Five sources are so confused by cloud material that identi- 
fying them as T Tauri disks or background sources is not 
possible with our data. The (lack of) variation and concentra- 
tion of the molecular emission rules them out as embedded 
YSOs, however. 

- All sources previously identified as disks using traditional 
methods such as Tboi and am are indeed recovered as disks as 
are all edge-on disk sources, identified by near-IR mapping 
(e.g. IRS 46, CRBR 2422.8-3423 and OphE-MM3). 

- Spectral line mapping of C'^O reveals that foreground lay- 
ers are predominantly present in the Ophiuchus ridge and 
the Oph-B core. Outside these regions, foreground layers are 
absent. These foreground layers are responsible for heavy 
reddening with typical Ay > 20, confusing classification 
schemes based on IR colors. 

The YSOs in the Ophiuchus clouds show a wide variety of 
characteristics, ranging from the rare Stage deeply embedded 
source VLA 1623 and Stage 1 embedded YSOs such as Elias 
29, IRS 43 and WL 12 to the embedded YSOs with wimpy 
envelopes in transition to the T Tauri stage like IRS 63. The 
results above imply that the molecular emission maps uniquely 
characterize the envelopes associated with embedded YSOs, 
while simultaneously identifying disks erroneously classified 
as embedded YSOs. Such characterization is needed before 
statistical studies can be used to estimate time scales and star 
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Fig. 17. Classification of sources using the spectral slope (upper left), bolometric temperature (upper right), HCO^ intensity (lower 
left) and HCO^ intensity combined with the bolometric temperature (lower right). In the upper two plots, embedded sources are 
shown with a diamond, disks are shown with a star, and confused sources with a triangle. In the lower two plots, diamonds are again 
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formation efficiencies within L 1688. 

One characteristic that has yet to be determined is the ratio of 
the disk to envelope mass. Interferometric studies with advanced 
facilities such as the SMA, CARMA or the PdB interferome- 
ter will allow the disks to be characterized and, through com- 
parison with continuum single dish studie s, the envelope (e.g. 
|j0rgensen et al.ll2005llLommen et al.ll2008l) . Such detailed infor- 
mation will be necessary as a starting point for the interpretation 
of future observations with Herschel and ALMA of the embed- 
ded YSOs in Ophiuchus and other clouds. A central role is envi- 
sioned for molecular line data. 
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Detailed source description 

GSS 30 and LFAM 1. These sources are both Stage 1 
sources with a separation of 15", but share a common cir- 
cumbinary envelope on scales up to 40". Strong HCO^ is seen, 
peaking in between both IR positions. Although outflowing 
gas was detected, the bulk of the emission is believed to be 
quiescent. 

VLA 1623. HCO+ is found to be extended on scales of 30", 
larger than that found for most sources. The lack of any emission 
at IR-wavelengths confirms this to be a Stage source. 
WL 12. This source is a good example of a small embedded 
Stage 1 YSO. It is embedded in a small ridge of material, 
extending from a north-east to west direction, as can be clearly 
seen in both dust and C'^O mapping. This ridge contributes 
only little to the HCO^ emission, which is concentrated at the 
source position. It is spatially resolved and peaks to 1 .4 K. This 
is much too brigth to be associated with a disk. 
LFAM 26. This Stage 1 is located at the tip of the Ophiuchus 
ridge that extends from LFAM 26 down to IRS 51. See § 5.2 
and appendix B for a discussion on the environment of LFAM 
26. The onset of this ridge is clearly seen, extending to the 
south-east direction, with much less material present in the 
north-west direction. Both HCO^ and 850 jum peak at the source 
position and are bright. Even if the emission originating in 
the cloud is subtracted from the HCO^ emission at the source 
position, it is too bright to be associated with a disk. This 
classification is confirmed by the very high of 1.27. 
WL 17. This is a Stage 1 source, but evolving towards a 
Stage 2, which we classified 'late' Stage 1. The emission of 
HCO^ is too bright to be associated with a disk. Since little 
to no extended emission is seen at scales of 40", the source is 
probably unobscured by cloud material. 

Elias 29. This is a well-st udied Stage 1 source dBoogert et alj 
I2OO2I iLommen et al." 2008). The analysis of HCO+ and contin- 
uum, both from SCUBA and SHARC-II, confirm this source to 
be embedded. 

IRS 37 and WL 3. The Stage 1 sources IRS 37 and WL 3, 

separated by 20" share a common envelope, with both C'^O and 
HCO^ peaking on both source positions. This envelope is highly 
non-spherical. The proximity of the two sources is probably 
responsible for this. Most of the mass is located near IRS 37, 
although the amount of HCO^ emission and concentration of 
the cores are equal. 

WL 6 and GY 256. The weak HCO+ seen for these two 
sources with a separation of 12" is brighter than that assumed 
to originate within a disk. In the dust maps, no emission is 
seen. WL 6 is probably a more evolved Stage 1 source with 
little envelope mass left. It is possible that GY 256 is a second 
component in a binary due to the rising spectra between 2 and 
24 jum. However, since the flux from GY 256 at 6 and 8 yum is 
less than at 4.5, it is more likely that GY 256 is a Stage 2 source. 
GY 256 is classified as confused. 

IRS 43. This is a Stage 1 source, with a well-defined envelope 
that is bright in HCO+ and C"^0. Both the SCUBA and HCO+ 
maps show a slightly elongated envelope with a possible second 
component about 30" to the east. This component is relatively 
much weaker in HCO^ than in the dust. However, no IR source 
is detected at that position. At the position of IRS 43, dense 
outflowing gas is prominently detected in HCO^ 4-3. 
IRS 44 and IRS 4 6. The Stage 2 source IRS 46 was stud- 
ied by iLahuis et al. I (I2OO6) in detail using IRS spectroscopy. 
Supplementary data from the JCMT confirmed that this source 
was a disk, reddened by a nearby envelope. The outflowing 



gas around the Stage 1 source IRS 44 is responsible for the 
reddening of IRS 46, as seen in Fig. [3] 

Elias 32 and 33. The Oph-B2 region dominates the environment 
around these two Stage 1 sources. At the position of Elias 33, 
HCO^ and 850 jum continuum clearly peaks. Elias 32 is located 
in the middle of the Oph-B2 region. A more pronounced 
SCUBA and a less pronounced HCO^ peak are seen within 10" 
of Elias 32 and aye. is in both cases near 0.0. 
IRS 54. Extended emission was found for SCUBA, HCO+ and 
C'^O. In all three cases, the emission is weak. IRS 54 is located 
outside of the main clouds. Classified as a late Stage 1, this 
source has accreted most of its envelope mass. 
IRAS 16285-2355. This is the only non-late stage 1 source that 
is located outside of the Ophiuchus ridge, Oph-B2 or Oph-A 
region. All emission in HCO^, C'^O and 850 jim continuum is 
associated with the source in a spherical envelope. 
IRS 63. This stage 1 source consists of a large disk with little 
envelope material. The observed HCO^ intensity of 1.2 K km 
s ' confirms that not all envelope material h as accreted onto the 
disk, as suggested bv lLommen et al. I (l2008h . 

CRBR 2324.1-1619. Due to the location of this source, 
within 2' of both VLA 1623 and LFAM 1/GSS 30, and located 
in between these bright sources, large amounts of material 
associated with the Oph A core are present, greatly confusing 
this source. However, the HCO^, C'^O or SCUBA emission do 
not peak at the source position. The spectrum is sharply rising 
with Q'iR=0.87. However, this is mainly due to the very high flux 
seen at 24 yum. 

GY 51. This source was not observed in HCO+. The SCUBA 
map shows no peak at the position, but bright extended emission 
in the direction of the Oph-A core. The source is classified 
as confused, but due to the ckir of 0.05, this source is likely 
a background disk, heavily extincted by the Oph-A cloud 
material. 

WL 19. In the continuum and molecular emission maps, a large 
amount of material is found, centered 20" north of WL 19. 
This material, probably associated with a pre-stellar core, is 
responsible for obscuring WL 19. Since the q-ir = -0.43 this 
source is not embedded. 

IRS 42. Most of the emission from IRS 42, in both 850 yum and 
HCO^ originates within the Ophiuchus ridge, near the core of 
IRS 43 and the location of CRBR 2422.8-3423.2. This source 
is classified as confused. However, with a q'ir=-0.03 it is more 
likely that this source is a disk, probably edge-on. The less 
likely possibility is that this source is a back-ground source. 

SSTc2d J162527.6-243648. This source, located far from 
the main L 1688 regions, is a disk. There is no detected SCUBA 
or HCO^ emission. The rising IR spectrum of 0.36, combined 
with the bolometric temperature of 105 1 K, is a strong indication 
that this source is an edge-on Stage 2 source. 

GSS 26. Extended SCUBA emission is seen near the source 
but is likely mostly associated with nearby cloud material. If this 
cloud emission is subtracted, a small unresolved core remains. 
The falling ffiR of -0.46 is the main reason that this source was 
classified as Stage 2. 

CRBR 2315.8-1700. Strong HCO+ emission and a strongly 
rising spectrum (q'ir=0.69) are seen for this source. The HCO^ 
map however suggests that this source is reddened by the large 
amount of material present in Oph A. An unresolved source 
is seen in the 850 jum data, but that core is not apparent in the 
HCO^ map. No variation was found in the HCO^ 4-3 at a level 
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Table 9. Source fluxes including IRAC 1-4, MIPS 1-2, SHARC-II, SCUBA and 1.3 mm data. 
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* 350/450 iim flux taken from lAndrews & William^ ilOOH) . Fluxes obtained with SHARC-II are given either in a 9" beam lAndrews & Wilhamsll2007h or in a 25" 
radius (This work)^ 

Data taken from lAndre etaj] <1993D . 

two sources within 4" with a total flux of 2e2 mJy 
" IRAC and MIPS-24 Fluxes marked with '*' are band-filled sources with a S/N<5 or sources with S/N<5. Such limits ai'e caused either by saturation or a very low 
emission. A 'K' added to the flux indicates that the obtained flux was flagged as 'Complex', with multiple detection within 2-3" 



of 0.2 K km s . Thus, this source is likely an disk, in line of 
sight of the cloud material. 

CRBR 2317.3-1925. CRBR 2317.3-1925 was observed with 
APEX-2a in HCO^ and was found to be reasonably strong with 
a peak temperature of 1.5 K. However, the source does not 
have a rising spectrum and the SCUBA core seen at the source 
position is very weak. In addition, the SCUBA emission is 
extended to the south-east. This source is classified as Stage 2. 
VSSG 1. Similar to CRBR 2317.3-1925, VSSG 1 does not have 
a rising spectrum. Unlike, CRBR 2317.3-1925, however, no 
HCO^ was found with APEX-2a. 850 fj.m emission is detected, 
but unresolved. The very high Tboi, caused by sharply falling 



spectra long-ward of 24 fim, confirms this source to be a Stage 
2 source, most likely edge-on. 

CRBR 2339.1-2032. This Stage 2 source has unresolved emis- 
sion at 850 fj.m. Although HCO^ was detected, the contribution 
is weak, only just above our cut-off for Stage 1, and does not 
seem associated with the source but with a patch of cloud 
material about 30" to the north-east of CRBR 2339.1-2032. 
WL 2. No HCO^ is detected associated with the source. In the 
dust, a small enhancement is seen at the source position, but it is 
not brighter than 0.15 Jy beam All HCO+ seen in the JCMT 
map is associated with a small core located between Oph-A and 
Oph-C. Combined with aiR = 0.02 and rboi=573 K, the source 
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Fig. 18. The individual IRS spectra (where available). Each spectra is plotted from 5 to 20 //m. All spectra are normalised to their 
own emission at 13.5 micron, indicated with a vertical dashed line. 



is classified as a disk, probably edge-on as well as extincted by 
the cloud material. 

GY 224. This source, located south of the Oph ridge, shows no 
rising spectrum and no HCO^. Since no 850 mum emission was 
found down to the noise limits of the COMPLETE map, this 
source is concluded to be a Stage 2. 

WL 20S. The emission seen in the SCUBA map does not peak 
at the position of WL 20S. In the SHARC-II map, a very weak 



(0.05 Jy beam"^) and unresolved peak is detected. In addition, 
no HCO^ was detected down to the limit adopted for disks. Its 
location on the Ophiuchus ridge confuses the environment of 
this source. The very high qtir = 2.7 suggests that this source 
must be an edge-on disk, even if located behind the Ophiuchus 
ridge. 

IRS 48. The disk source IRS 48 is isolated, away from the 
main cores in Ophiuchus. A weak, unresolved peak is seen in 
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the SCUBA. The absence of any C"*0 and HCO+ confirms the 
classification to be Stage 2. 

GY 312. GY 312 is not detected in the dust or gas mapping. In 
additon, no contribution from large-scale material is seen. Its 
clearly rising IR spectrum identifies this source as an edge-on 
Stage 2 source. 

IRS 51. IRS 51 is located at the south end of the Ophiuchus 
ridge. The ridge is clearly seen in SCUBA as well as in C"^0 
and HCO^ maps around IRS 5 1 . Although a clear peak is seen 
at IRS 5 1 in the dust map, no peak was seen in the molecular 
line mapping at the position of IRS 51. It is concluded that IRS 
5 1 is a Stage 2 source with strong disk emission, reddened by 
the Ophiuchus ridge. 

SSTc2d J162741.6-244645. This stage 2 source, located south 
of the Ophiuchus ridge, was not detected in either HCO^ or 
SCUBA. The rising spectrum suggests that this source must be 
edge-on since no cloud material is nearby to redden the source. 
SSTc2d J162748.2-244225. For this source, located east of 
the southern tip of the Ophiuchus ridge, no SCUBA or HCO^ 
emission was found. The very high cir of 1 .98 suggests that this 
source must be an edge-on disk. 

SSTc2d J162857.9-244055. This source is located between the 
L 1688 and L 1689 regions. The combination of a steeply rising 
spectrum and an absence of HCO^ and 850 jim continuum 
indicates that this source is a edge-on Stage 2 source. 
Haro 1-4. Located north of the main Ophiuchus clouds, this 
source has no HCO^ or 850 //m associated with it and clearly is 
a disk. 

DoAr 25. Although a bright peak is seen in the 850 jum map, 
the (TiR of -1.12 clearly confirms this to be a disk source. 
OphE MM3. This source, located within t he Ophiuchus ridge , 
was first discovered in the mm-regime by iMotteetalJ (Il998l) . 
Both the 850 /im continuum and HCO^ show little material 
associated with an env elope at this source. Near-IR observations 
jBrandner et al.ll2000l) confirm it to be an edge-on disk. 
SR 21. This source has a very high bolometric temperature of 
1070 K. Although a clear signal is seen in the 850 jim data, no 
emission is seen in HCO ^. It has been confirmed as a cold disk 
with a large inner hole bv lBrown et alj ( |2007[) . 
CRBR 2422.8-3423. CRBR 2422.8-342 3 is located along the 
Ophiuchus ridge. From near-IR imaging (jBrandner et al.ll2000l 
iPontoppidan et alJ l2005h it is known that CRBR 2422.8-3423 
is an edge-on disk and located in line of sight of the IRS 
43 envelope. Even at far-IR wavelengths, emission from the 
Ophiuchus ridge rather than the source dominates. HCO^ 
mapping confirms this. 

SR 9. This source is not detected at 850 jum. Combined with 
^iR < -1, it is concluded that SR 9 is a Stage 2 source. 
2MASS 16282. Although 2MASS 16282 has a very low 
bolometric temperature of 66 K, the q-ir < -1 identifies this 
source as a disk. 



